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Abstract 

We consider a model of quantum-wire junctions where the latter are described by conformal- 
invariant boundary conditions of the simplest type in the multicomponent compactified mass¬ 
less scalar free field theory representing the bosonized Luttinger liquids in the bulk of wires. 
The boundary conditions result in the scattering of charges across the junction with non¬ 
trivial reflection and transmission amplitudes. The equilibrium state of such a system, cor¬ 
responding to inverse temperature /3 and electric potential V, is explicitly constructed both 
for finite and for semi-infinite wires. In the latter case, a stationary nonequilibrium state 
describing the wires kept at different temperatures and potentials may be also constructed 
following Ref. [32]. The main result of the present paper is the calculation of the full counting 
statistics (FCS) of the charge and energy transfers through the junction in a nonequilibrium 
situation. Explicit expressions are worked out for the generating function of ECS and its 
large-deviations asymptotics. For the purely transmitting case they coincide with those ob¬ 
tained in Refs. [10, 11[, but numerous cases of junctions with transmission and reflection are 
also covered. The large deviations rate function of FCS for charge and energy transfers is 
shown to satisfy the fluctuation relations of Refs. [2, 12[. The expressions for FCS obtained 
here are compared with the Levitov-Lesovic formulae of Refs. [29, 28[. 


1 Introduction 

The transport phenomena in quantum wires (carbon nanotubes, semiconducting, metallic and 
molecular nanowires, quantum Hall edges) and, in particular, across their junctions, have at¬ 
tracted a lot of interest in recent times, see e.g. [16, 14]. To a good approximation, the charge 
carriers inside the wires may be described by the Tomonaga-Luttinger model [46, 42, 22, 26]. In 
the low energy limit, such a model reduces to a relativistic 1+1 dimensional interacting fermionic 
field theory that can also be represented by free massless bosonic fields. The junction between 
the leads couples together the conformal field theories (CFTs) describing at low energies the bulk 
volumes of the wires. Specific features of the coupling depend on how the junction is realized. 
Various models that couple two or more wires locally at their connected extremities were consid¬ 
ered in the literature, see e.g. [20, 36, 34, 35] where important results about transport properties 
of such models of wire-junctions were obtained. The low-energy long-distance effect of the inter¬ 
action at the junction may be described with the use of boundary CFT, similarly as the effect of a 
magnetic impurity in the multi-channel Kondo problem [1]. Even if the coupling of the Luttinger 
liquid theories introduced by the junction breaks the conformal symmetry, the latter should be 
restored in the long-distance scaling limit. In the scaling limit, the effect of the junction will be 
represented, using the “folding trick” of ref. [49], by a conformal boundary defect in the tensor 
product of the bulk CFTs of individual wires [7]. Such a boundary defect preserves half of the 
conformal symmetry of the bulk theory. Examples of conformal boundary defects that describe 
the renormalization group fixed points of Luttinger liquid theories with a coupling localized at 
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the junction were discussed in [20, 36, 34, 35]. It was also realized that the boundary CFT de¬ 
scription of the junction of wires gives via the Green-Kubo formalism a direct access to the low 
temperature electric conductance of junctions [35, 39, 40] that measure small currents induced 
by placing different wires in slightly different external electric potentials. Getting hold of the 
transport properties of the quantum-wire junctions beyond the linear response regime is more 
complicated, see [20] for an early result using an exact integrability of a model of contact between 
two wires. The GFT approach seems also helpful here. It was shown in [10, 11, 12, 19] that for 
some boundary defects (those with pure transmission of charge or energy), not only the electric 
and thermal conductance but also the long-time asymptotics of the full counting statistics (FGS) 
of charge and energy transfers through the junction may be calculated for the wires initially 
equilibrated at different temperatures and different potentials. Moreover, steady nonequilibrium 
states obtained at long times from such initial conditions could be explicitly constructed. Physical 
restrictions for the applicability of the GFT approach in such a nonequilibrium situations were 
also discussed in some detail in those works, in particular in [11], see also [8, 18, 4, 15]. The 
incorporation of junctions corresponding to boundary defects with transmission and reflection 
into that approach poses more problems, although for a junction of two GFTs a general scheme 
has been recently laid down in [13], together with some examples. 

The present paper arose from an attempt to calculate the FGS for nonequilibrium charge 
and energy transfers for simple conformal boundary defects with transmission and reflection. 
We describe each of N wires by a compactified free massless 1 -|- 1-dimensional bosonic field, 
with the compactihcation radius related to the Luttinger model coupling constants that may be 
different for different wires. The product theory is a toroidal compactihcation of the massless 
^component free held, i.e., on the classical level, its held takes values in the torus 1/(1)^. In 
such a theory, we consider the simplest conformal boundary defects that restrict the boundary 
values of the held at the junction to a subgroup B C U{1)^ isomorphic to the torus U{\)^ 
with M < TV. In the string-theory jargon, B is called the D(irichlet)-brane [38]. First, we study 
the wires of hnite length L with the rehecting boundary condition at their ends not connected 
to the junction. The overall [/(l)-symmetry of the theory is imposed, leading to the conservation 
of the total electric charge. We show that the boundary defect gives rise to an TV x TV scattering 
matrix S that relates linearly the left-moving and the right-moving components of the electric 
currents in various wires. The classical theory described above may be canonically quantized 
preserving the latter property. The exact solution for the quantum theory includes the formula 
for the partition function of the equilibrium state corresponding to inverse temperature /3 and 
electric potential V and for the equilibrium correlation functions of the chiral components of the 
electric currents. The thermodynamic limit T —^ oo may then be performed giving rise to a 
free-field theory that was constructed directly for L = oo in [32]. In that limit, the equilibrium 
correlation functions involving only left-moving (or only right-moving) currents factorize into the 
product of contributions from the individual wires. This property was used in [32], following the 
earlier work [31], to construct a nonequilibrium stationary state (NESS) where the correlation 
functions of left-moving currents factorize into the product of equilibrium contributions from 
individual wires, each corresponding to a different temperature and a different potential. The 
NESS correlation functions involving also the right-moving currents are reduced to those of the 
left-moving ones using the scattering relation between the chiral current components. Following 
the approach of [10, 11], we show that such a state is obtained if one prepares disconnected wires 
each in the equilibrium state at different temperature and potential and then one connects the 
wires instantaneously and lets the initial state evolve for a long time [41]. 

The main aim of the present paper is the study of the FGS for charge and energy (heat) 
transfers through the junction modeled by the brane defect of the type described above. Similarly 
as in [11], the FGS is obtained from a two-time measurement protocol. First, the total charge and 
total energy is measured in each of the disconnected wires of finite length L prepared in equilibria 
with different temperatures and potentials. Next the wires are instantaneously connected and 
evolve for time t with the dynamics described by the held theory with the brane defect. After 
time t, the wires are disconnected again and the second measurement of total charge and total 
energy in individual wires is performed. The FGS is encoded in the characteristic function of the 
probability distribution of the changes of total charge and total energy of individual wires. The 
above protocol is not practical for long wires as the total charge and and total energy of the wires, 
unlike their change in time, behave extensively with L, but a similar charge and energy transfer 
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statistics should be obtainable from an indirect measurement protocol where one observes the 
evolution of gauges coupled appropriately to the wires and registering the flow of charge and 
energy through the junction, see [29, 30]. In our model, we compute the generating function of 
FCS of charge transfers explicitly for any L and t and confirm that it takes for large t the large- 
deviations exponential form that is independent of whether L is sent to infinity first or, e.g., kept 
equal to t/2. The equality of the large deviation forms for the two limiting procedures appears, 
however, to be less obvious than one could have expected. The choice L = t/2 leads to the 
simplest calculation of the large deviation rate function and was implicitly employed in [10, 11], 
where it was argued that it reproduces correctly the large deviations of the FCS for the junction 
of semi-infinite wires. We also compute explicitly the generating function of the FCS for heat 
transfers for L = t/2 and its large deviations form. The case of general L and t could be also 
dealt with but the corresponding formulae are considerably heavier and we did not present them 
here. The generating function of the joint FCS of the charge and energy transfers for L = t/2 and 
its large deviations form were also obtained. To our knowledge, the calculations of FCS presented 
in this paper are the first ones obtained for junctions with transmission and reflection modeled by 
conformal boundary defects. It should be mentioned, however, that in a different physical setup, 
the FCS of charge transfers across an inhomogeneous Luttinger liquid conductor connected to 
two leads with distinct energy distributions was obtained by a “nonequilibrium bosonization” in 
[24, 25, 33]. 

The present paper is organized as follows. In Sec. 2, we briefly recall the description of rela¬ 
tivistic free massless fermions and bosons on an interval. We discuss the correspondence between 
the two theories and how it extends to the case of the Luttinger model of interacting fermions. 
Sec. 3 describes in detail the model of a junction based on a toroidal compactification of the 
multi-component massless bosonic free field with a boundary defect of the type mentioned above. 
We discuss first the classical theory on a space-interval of length L and subsequently canonically 
quantize that theory in Sec. 4. In particular, we show how the scattering matrix S relating the 
chiral components of the electric current arises from the brane describing the boundary defect. 
Sec. 5 constructs the equilibrium states of the quantized theory labeled by inverse temperature 
/3 and electric potential V. In Sec. 6, we discuss the Euclidean functional integral representation 
of the equilibrium state and in Sec. 7, its dual closed-string representation resulting from the 
interchange of time and space in the functional integral. The closed-string picture is particularly 
convenient in the thermodynamic limit T —>■ oo of the equilibrium state that is analyzed in Sec. 8. 
Sec. 9 discusses the NESS of the junction of semi-infinite wires kept in different temperatures and 
different electric potentials. By considering the nonequilibrium state for close temperatures and 
potentials, we obtain as a byproduct the formulae for the electric and thermal conductance of 
the junction. The central Sec. 10 is devoted to the analysis of FCS for charge and heat transfers 
through the junction. Subsecs. 10.1 and 10.2 treat the charge transport, Subsec. 10.3 that of heat, 
and Subsec. 10.4 the joint FCS for both. Sec. 11 compares the generating function of FCS for 
charge and heat transfers obtained in this paper with those given by the Levitov-Lesovik formu¬ 
lae for free fermions [29, 30] and free bosons [28]. In Sec. 12, we specify our general formulae to 
few simplest cases of junctions of two and three wires. Finally, Sec. 13 collects our conclusions and 
discusses the possible generalizations and open problems. Appendix A contains the calculations 
of the generating functional of FCS for charge transfers at general t and L. Appendix B per¬ 
forms the computation of certain bosonic Fock space expectations that are needed to obtain the 
generating function of FCS for heat transfers through the junction. Appendix C calculates the 
quadratic contribution to the Levitov-Lesovik large-deviations rate function of charge transfers 
for free fermions. 

Acknowledgements: The authors thank D. Bernard for discussions on nonequilibrium CFT 
and J. Germoni for Ref. [44]. A part of the work of K.G. was done within the STOSYMAP project 
ANR-ll-BSOl-015-02. 

2 Field theory description of quantum wires 

2.1 Classical fermions 

Consider a fermionic 14-1-dimensional field theory describing noninteracting conduction electrons 
in a quantum wire of length L. To a good approximation such electrons have a linear dispersion 
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relation around the Fermi surface. For simplicity, we shall ignore here the electron spin. The 
classical action functional of the anticommuting Fermi fields of such a theory has the form 


L 

^ J dt J dx, 

0 


where d± = ± dx), with the boundary conditions 






( 2 . 1 ) 


( 2 . 2 ) 


We use the Fermi velocity vf to express time in the same units as length. The classical equations 
obtained by extremizing action (2.1) are 

= 0 = , a+V'’' = 0 = a+vi’' (2.3) 

and their solutions take the form: 

^ CpC = V’’'(i,-a;), (2.4) 

^ (2.5) 


The space of classical solutions comes equipped with the odd symplectic form 

L 

A X.l/ I A ^ A 

L.p 


fl = ^ / [S-ip^ A 6ip^ + S'ip'^ A Stp'^'\ dx = i 6cp A Sc 
0 p^^+i- 


leading to the odd Poisson brackets 


{cp,Cp^} — 0 — {cp,Cp'} , {cpjCp'} — i (5p 


The U{1) symmetry 


corresponds to the Noether current 
70 _ 1 tj.i 


^ gia^Ar- 




with the chiral components 


/ = i(jo - Jl) = , r = UJ^ + C) = 


1 7 tO , 7 l\ _ 1 


and the conserved charge 


The classical Hamiltonian is 


clc„. 


Q = y j°(t,x)dx= L-pL-p 

0 


H = 


[il^^dxip^ - ^p^^x'^p"'] dx = ^ ^ k CpCp . 

pC^+-^ 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

(2.12) 


2.2 Quantum fermions 

Quantized Fermi helds and 'ip'' are given by expressions (2.4) with operators Cp and their 
adjoints c|, satisfying the canonical anticommutation relations 

[cpjCp']^ = 0 = [Cp,Cp,]^, [cp,Cp,]^ = Sp^pi . (2.13) 
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They act in the fermionic Fock space Jy built upon the normalized vacuum state |0)j: annihilated 
by Cp and cLp for p > 0 (the annihilation operators of electrons and holes, respectively). Upon 
quantization, fields and •0’’ become the hermitian adjoints of and ip'". The quantum 
C/(l) currents have the chiral components 

/ = r = ^:rr-- (2.14) 

and the conserved 17(1) (electric) charge is^ 

L L 

Q= f J°dx= f{J^+ r)dx = ■ 4^p ■ ■ ( 2 - 15 ) 

0 0 P6Z+T 

The fermionic Wick ordering putting (electron and hole) creation operators Cp and for 
p < 0 to the left of annihilators Cp and cLp for p > 0, with a minus sign whenever a pair is 
interchanged, assures that the vacuum |0)j. has zero charge. The quantum Hamiltonian is 

^ = t( 12 P ■ 4^p ■ - (2.16) 

where the constant contribution is that of the zeta-function regularized zero-point energy 

oo n 

+ |C(-1) =-4- (2-17) 

p<.0 n—1 n—1 

2.3 Classical bosons 

Consider now a bosonic 1+1-dimensional massless free field if{t,x) defined modulo 27r on the 
spacetime M x [0, L], with the action functional 

L 

-S'M = 4^ j dt j dx. (2.18) 

0 

We shall impose on ippt, x) the Neumann boundary conditions 

dMt,0) = 0 = dMt,L)- ( 2 . 19 ) 

Such a scalar field will be viewed as having the range of its values compactified to the circle of 
radius r with metric The classical solutions extremizing action (2.18) have the form 

(p{t,x) = (p^{t,x) + ip^{t,x) (2.20) 

with 

4{t, x) = ipo + ^ao(l + a;) + i ^ 5 ^ '^ 2 n e ^ = p’'(f, -x) (2.21) 

and d 2 n = C(- 2 n- The labeling of modes a 2 n by even integers is for the later convenience. The 
symplectic form on the space of classical solutions is equal to 

n = ^6ao A Sipo - ^ A 6a-2n , (2.22) 

n^O 


leading to the Poisson brackets 


+o} — 2r , {rr2n,C^2n'} — ‘2tI\T Sn-\-n'^0- 


(2.23) 


The [/(I) symmetry 


ip 1-^ (fi + a 


(2.24) 


'^Here and below, we measure the electric charge in the negative units — e so that electron’s charge is -|-1. 
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corresponds to the Noether current 


with the chiral components 

= ^d±(p{t,x) = 


^2 > _ 7rin(t±x) 

iT “2n e I , 




( 2 . 25 ) 


( 2 . 26 ) 


where the upper sign pertains to the left-moving component depending on x~^ = t + x and the 
lower one to the right-moving one depending on a;“ = t — x. The classical Hamiltonian takes the 
form 

L 

^ ^ + {d^ipf]dx = ^^a2na-2n- ( 2 . 27 ) 

n riGZ 


2.4 Quantum bosons 

The space "Ho of quantum states corresponding to the zero modes may be represented 

as L^([/(l)), with viewed as the angle in U{1). ao acts then as —2ir“^g^ assuring the 
commutation relation [ao, po] = —2ir“^. An orthonormal basis of "Hq is composed of the states 

\k) = (2.28) 

with 

ao|fc) = 2r-2fc|fc). (2.29) 

The excited modes a 2 n = ct^- 2 n with the commutation relations 

[o2n, 02n'] — 2tIT 6n-i^n' ,0 (2.30) 

are represented in the bosonic Fock space built upon the vacuum state |0)t, annihilated by 
Q! 2 n with n > 0. The total bosonic space of states is "Hf, = "Ho 0 J^b- We shall identify "Ho with 
its subspace "Ho^lO);, and the state |0) G "Ho with the vacuum |0)(8)|0)b. The chiral components 
J^’’’ of the quantum 17(1) current are given by the right hand side of Eq. (2.26). The conserved 
U(l) charge takes the form 

L 

Q = y J° dx = |x^ao (2-31) 

0 

so that Q\k) = k\k) and it acts trivially in J^b- The quantum Hamiltonian requires a bosonic 
Wick reordering putting the creators a_ 2 n for n > 0 to the left of annihilators a 2 n in the classical 
expression. Explicitly, 

00 

» = + P-32) 

n—1 

where the constant term is the contribution of the zeta-function regularized zero-point energy 

00 

^E2^"'' = ^C(-1) = -^. (2.33) 

n—1 

2.5 Boson-fermion correspondence 

In one space dimension there is an equivalence between quantum relativistic free fermions and 
free bosons that provides a powerful tool for the analysis of such systems [43, 27], see also [45] 
for the historical account. In the context of the fermionic system described in Sec. 2.2, such an 
equivalence involves the free bosonic field of Sec. 2.4 with the compactification radius r = and 
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is realized by a unitary isomorphism X : 'H.b —>■ "H/ that maps vacuum to vacuum, I|0)& = |0)/, 
and intertwines the action of C/(l) currents and the Hamiltonians^. In particular, 

I otiu = ^ : clcp+n ■ I. (2.34) 

The Fermi fields are intertwined by X with the bosonic vertex operators: 

= I : 


= I 


:(*+^) p-ivo p-^aoit+x) 


-^“2716 

X e"<“ 


7rin(t + a3) 

- Z - 


■^“2ne 

0 n >0 


n(t + cc) 


i}\t,x)X = X : e^^v\t,x) 


X 



^ivo 


g^aod+x) 


“ Yl w“2ne 

X e 


7Tin(t + x) 

Z 


- 12 i“2ne 
e 


7rin(t + £c) 

—^ 


(2.35) 


2.6 Luttinger model 

The interaction of electrons near the Fermi surface gives rise to the addition of a perturbation 
to the free field Hamiltonian (2.16) that in the leading order takes the form of a combination of 
quartic terms in the free fermionic fields: 

L 

[ 2 ff 2 (: y/ :)(: 0 + 54 ((: O" + (: VV O")] dx + const., (2.36) 

0 

where an infinite constant is needed to make the operator well defined in the fermionic Fock 
space. Such a perturbation defines the Luttinger model of spinless electrons in one-dimensional 
crystal [46]. The crucial fact that enables an exact solution of such a model is that, under the 
bosonization map, the above perturbation becomes quadratic in the free bosonic field: 

L 

X=X ^ J : [ 2 g 2 {d+(p){d_(p) + 54 iid+(pf + {d^ipf)] : dx -f const. 

0 

L 

= ^ 4^ / • [(54 + g2){dtipf + (54 - g2){dx^f] ■ dx -f const., ( 2 . 37 ) 

0 

where on the bosonic side the Wick ordering takes care of the diverging part of the constant on 
the fermionic side. The perturbed bosonic Hamiltonian has then the form 

L 

= H-\- ^ J : [(27r-\-g4-\-g2){dt(p)‘^-\-{27r-\-g4-g2){da:(f)‘^] : dx + const. (2.38) 
0 

in terms of the free field (p(t, x) with the compactification radius r = ^/2. corresponds to 

the classical Hamiltonian 

L 

= i y [(2^ + 34 + 92)11^ + [dxvf] dx, (2.39) 

0 

where n(t,x) = ^{dtip){t,x) is the held canonically conjugate to (p{t,x). The classical La- 
grangian related to the above classical Hamiltonian is obtained by the Legendre transform: 

L a-^L 

= i/ d* = £ / [ia,~pf - at' (2.40) 

0 0 

^We choose to represent free fermions by bosons compactified on the radius r = \/2 rather than on the more 
frequently used dual radius r = as better suited to the fermionic boundary conditions (2.2). 
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for X = ax', where 


fL = T{ — / ^'^+ 94-92 

2 - “ Y 2ii-+g4+S2 ’ 


^ _ 'h’ren 

~ VF 


\/ (27r+g4)^-si 
Sir 


(2.41) 


Hence after the change of the spatial variable, Lagrangian becomes that of the free bosonic 
field compactihed on the radius r that is different from r = \/2 if g 2 ^ 0. The factor a gives 
the multiplicative renormalization of the wave velocity vp due to the interactions (we assume 
that |27r + 34 ! > |(/ 2 |)- The quantization of the free bosonic theory compactihed at radius r 
discussed in Sec. 2.4 provides the exact solution of the Luttinger model on the quantum level. 


3 Bosonic model of a junction of quantum wires 


In the spirit of the “folding trick” of [49, 37], see Fig. 1, we shall model a junction of N quantum 
wires by a compactihed free held g(t,x) with A^-components gi(t,x) = € U(l) dehned 

on the spacetime K x [ 0 , L], with the action functional 



) dx 


(3.1) 


and appropriate boundary conditions. The compactihcation radii may be different for different 
wires, corresponding to different quartic coupling constants g 2 i and ga in the Luttinger models 
describing the electrons in the individual wires, see Eq. (2.41) of Sec. 2.6. We shall impose the 
Neumann rehecting boundary conditions at the free ends of the wires: 


a,<^(f,L)=0, (3.2) 

where ip = {pi). Note that we use the rescaled spatial variables in the wires so that the lengths 
of the wires in physical variables are hxed to UiL. This will not matter much because the length 
L will be ultimately sent to inhnity. 



0 L 


boundary 

defect 


Folding trick 


The “boundary defect” representing in the folding trick the junction of wires at x = 0 will 
be described by the boundary condition requiring that the 17(1)^-valued field g belongs to a 
^ r ^ * 

g{t,G)€B = n{U{l)^)cU{l)^, (3.3) 

where k : 1/(1)'^—> 1/(1)'^ is a group homomorphism 

M _ N 

^ (3-4) 


specified by integers k™. We shall assume that k is injective so that k{U{1)^) = U{1)^. As 
may be seen from the Smith normal form of matrix (k™), such a property is assured if and only 
if the M X N matrix (k™) has rank M and the g.c.d. of its M x M minors is equal to 1, see 
Proposition 4.3 of [44]. In particular, M < N necessarily. Consider matrices T = (T'”^'" ) and 
P = {Pii') defined by the relations 


N 


M 


rpmm' 2 m m' P 






(3.5) 









Matrix P defines the projector on the snbspace of spanned by the vectors = (k™, •.. 
that is orthogonal with respect to the scalar product 

a-b = '^rfaibi (3.6) 

i 

in . The boundary condition (3.3) implies that 

P^dMt,0) = 0, (3.7) 

where P-^ = I — P. The stationary points of the action functional (3.1) satisfy, besides the 
imposed boundary conditions, the equations 

i9t = 0, (3.8) 

PdMt,0) = Q- (3.9) 


Note that relations (3.7) and (3.9) imply mixed Dirichlet-Neumann boundary conditions at a: = 0 
for massless free fields The solutions of the classical equations decompose in terms of 

the left- and right-movers: 


x) = ip^{t -I- a;) -I- (p^{t - x) 


(3.10) 


with 






2L ^0 


(tiLx)+i 


E: 

n^O 


^(t±:x) 


where the upper sign relates to the and the lower one to tp’’, and 


(3.11) 


a 




— OCr 


= {p- p^)a: 


Pa 


i,r 

2n+l 


= 0 , 


P^ai^: = 0 


(3.12) 


with real ipo = such that e“^“ G B. In particular, al^^+i = —(^ 2 n+i ^ 2 n = ^ 2 n 

The space of classical solutions comes equipped with the symplectic form 


fl = ^ai ■ ^ (3.13) 

n^O 


which determines the Poisson brackets of functionals on that space that may be directly quantized. 

The particular case when n in (3.4) is the identity mapping of U{1)^, corresponding to the 
“space-filling” brane Bq = U{1)^, describes the disconnected wires. In this case, P = I (i.e. P 
is the identity matrix) and field (p satisfies the Neumann boundary conditions both at x = 0 
and X = L and only the even modes = *^ 2 n appear. One obtains in this case the product 
of N theories considered in Sec. 2.3. 


4 Quantization 

4.1 Space of states 

The quantization of the bosonic theory of Sec. 3 is again straightforward but a little more involved 
than for the disconnected wires. Let us first quantize the zero modes. According to the boundary 
conditions, 

«o = «o = <Po = 

m m 

where {tpm), to = 1 ,...,M, are angles parameterizing 17(1)'^, so that 

i ^ A JV'm' • (4.2) 

m,m' 

The corresponding Poisson brackets are 

{/3m, V'm'} = -2(T-i)w (4.3) 
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leading to the commutators 


[(3^, V'm'] = ■ (4.4) 

Keeping in mind that the angular variables tprn are multivalued, the above commutators will 
be represented in the Hilbert space Ho = of functions of M angles 'ipm, square 

integrable in the Haar measure, by setting 

/3m = ~2i^^(T )mm' ' (4-5) 

m' 

An orthonormal basis of "Ho is given by the states 


such that 


= 2 

? 

For the excited modes, it is convenient to i 
in such that 

^3 ■ ~ ’ ^-^3 ~ -^3 

and the projected modes 

c 

with the inverse formulae 


1 

(4.5) 

for fc™ G Z 

(4.6) 


(4.7) 

\k\..k^). 

(4.8) 

Vj)^i of vectors A^ = (A^i,. 

.., ^jn) 

= Kj for j > M 

(4.9) 




N 

i=i 


(4.10) 

(4.11) 


Note that relations (3.12) imply that a( 2 n+i) j ~ J < Af and < 5 ^ 2 n) j 4 ^ 

Poisson brackets of the non-zero operators oinj take the form 




leading to the commutators 


— n6jjiSn+n',0 ■ 

In the standard Fock space quantization, we take 


M 

3F e — ® 3F Qj , 
1=1 


(4.12) 

(4.13) 

(4.14) 


1=1 j=M+l 

where Hej and Hoj are generated by vectors 

«-(2«i)j ■••ai(2 ni)j |0)e ^^0, ni>...>ni>l, j e M}, (4.15) 

di(2„^+i)j ... di(2„,+i)j |0)„ / > 0, ni > ... > n; > 0, j G {M -f 1,..., iV} (4.16) 

with the scalar product determined by the relations 

4n)il0)e = 0 for n>0, 4„+i),|0)o = 0 fo^ «>0, (4.17) 

e(0|0)e = o(0|0)o = 1 , (dfo)^ = . (4.18) 

The total Hilbert space of states of the theory is 

TL = Ho 0 0 Ho (4.19) 

and in the following we identify 

\kH..k^)=\kH.. k^) ® |0)^ 0 |0)„ (4.20) 
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4.2 Currents, charge and energy 

We shall be interested in the system that possesses global U{1) symmetry acting on fields by 
{gi{t,x)) I—>■ {ugi(t,x)) for u G U{1). Invariance of the theory requires that this action preserves 
the brane B = k{U{M)) C U{N). This holds if and only if the vector 1 = (1,..., 1) is in the 
image of projector P, i.e. if PI = 1 or 

= 1 for fol fo (4-21) 

i' 

The Noether (electric) current corresponding to the U(l) symmetry has then the form 


N 


N 




in terms of the currents in individnal wires with the left-right moving components 

_ T-Ctlbx) 

^ni ^ 


= \{J° ^ J^){t,x) = ^\(dt±d^)ip,{t,x) = 


(4.22) 


(4.23) 


defining jf{t,x) and J^{t,x) as functions, respectively, of < + x and t — x for any real t and x. 
We shall use formulae (4.23) also for quantum currents. At x = 0, the left and right currents are 
linearly related: 

Jl(t,Q) = Y^S,e4{t,Q), (4.24) 

where 

= P^'^ - P^h (4.25) 

according to (3.12) and the explicit expression (3.5) for the matrix P. The N x N ”S'-matrix” 
S = {Sii' ) describes the flow of the currents through the junction of wires. It satisfies the relations 


Si'I = r, 



^ Sw = 1 , 


In other words, 

Sii' — XiOii'Tj^, 

where O = {Ow) is a symmetric orthogonal matrix such that 


'^TiOii' = n, . 
i 


(4.26) 


(4.27) 


(4.28) 


We shall use matrices S and O interchangeably. For N = 2, there are two possibilities: 


O = 


il ?) 


or 


O = 


(1 


2rir2 


2 X1X2 
„2 „2 
"2 ~ X 


?)■ 


(4.29) 


The first case corresponds to the identity embedding k describing the disconnected wires whereas 
the second one corresponds to the diagonal embedding of U{1) into U{\Y that leads to a 
nontrivial junction. In the last case, the xi = r 2 case corresponds to off-diagonal matrix O = S 
with unit non-zero entries, i.e. to the pure transmission of currents through the junction, but for 
xi X 2 the currents are partly transmitted and partly reflected at the junction. 

Eq. (4.24) implies that the right currents are linear combinations of left currents if considered 
as functions of real t and x: 


4 (t: 4 =Y1 ■ 

i' 

At X = L, i.e. at the ends of the wires, the left and right currents are equal: 

4it,L) = 4{t,L) 


(4.30) 


(4.31) 
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which implies that 


Jl{t,x) = Jl{t,-x + 2L) (4.32) 

if we treat the currents as functions of real t and x. The quantum currents satisfy the equal-time 
commutation relations 

[Jf(t,x), 4{t,y)] = 51 E - P^4) ^'4 -y + 2nL) 

= -[J4t,x),J4t,y)], (4.33) 

[4it, x), 4 (t, J/)] = 5J E - ^«')) S'ix + y + 2nL), (4.34) 

In particular, the left-moving currents commute among themselves at equal times if their positions 
do not coincide modulo 2L. Similarly for the right-moving currents. The left-moving currents 
commute with the right-moving ones at equal times if their positions are not opposite modulo 
2L. Note that for 0 < x,y < L the only terms that contribute to (4.33) and (4.34) have n = 0 or 
n = 1, respectively, so that for such values of x and y the commutation relations of currents do 
not depend on the choice of brane B. This permits to identify for different junctions the algebras 
of observables generated by currents 4’^{0,x) with 0 < a; < T, i.e. localized away from the 
contact point. In particular, we may identify such observables for disconnected wires with those 
for connected wires, with the physical meaning that their measurement just before and just after 
establishing or breaking the connection between the wires should give the same result. Whatever 
the junction, the total charge 

N 

Q{t) = Y,Q4), (4.35) 

where Qi{t) are the charges in the individual wires, 

L 

Qi{t) = j J°(^,x)dx, (4.36) 

0 

is conserved: 

N ^ 

= E /4 + dt4 a;))dx 

i=i I 

N ^ N 

= E / 4 - 9x44 x)4x = E [44 4 - 44 a:)] 

i—1 Q i—1 

N 

= 2 E - P^'^)4' 4 0) = 0 (4.37) 

2,z' = l 

due to (4.21). In terms of the modes, 

e = E I• (4.38) 

i 

Operator Q acts only on "Ho : 

=E E \4 .. .k^) = {p,T~^k)\k4 . .k^) , (4.39) 

i m,m' 

where p = (p™) G with 

p^ = Y^4nT = l-K^, (4.40) 

i 

and 

(a,b) = E«'”^’”> (4-41) 

m 
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denotes the standard scalar product on , to be distinguished from the one of (3.6) used in 
Note that the spectrum of Q is composed of integers, as must be the case for the generator 
of a unitary action of U{1) group. Indeed, since for each 1 < z < iV, 

M 

Y. = {P^h = 1 ( 4 . 42 ) 


by (4.21), the injectivity of the homomorphism (3.4) implies that the sums ^p"^{T ^)mm' are 

m 

integers. This is not the case for (non-conserved) charges in the individual wires 


L 

Qi{t) = J dx 

0 



irf 1 £ -7ri(2n + l)t 

~ Z_^ 2n+l ‘^(2n+l)i ® ^ 

n 


(4.43) 


The energy of the bosonic system of Sec. 3 is given by its classical Hamiltonian that may be 
expressed in terms of the left and right moving currents by the formula 


N 

= [ (('^/(^’^))^ + ('^r(^>a;))2) dx. 


(4.44) 


Its conservation, that holds independently of the condition (4.21), results from the identity 


® = E 27 rr-^ [{ 4 nt,x) - {Jir{t,x)Y 


dt 

whose right hand side vanishes because 

{jrfit,L) = {jfr{t,L) 


=L 

0 


(4.45) 


(4.46) 


for all 1 < i < A^, and because 

N N N 

Y ^■"^»'^«”4(t,0)4.(t,0)=^r-2(j/)2(t,0) (4.47) 

i—1 2,2',^”=! 2=1 

in virtue of (4.26). The quantum Hamiltonian H is given by the Wick reordered version of the 
classical expression: 


N 

P = Y^ /(■ xf ■■ + ■ X {t, xf ■■)dx + I 


TT N—3M 


48 


M 


= 4l «o ■ "0 + 5r “(- 


2n)i^{2n)i 


2=1 n>0 


N 

-£ ~£ , N—3M (A AQ\ 

2L '^(-(2n+l))i'^(2n+l)i 48 ’ (4-4§) 

n>0 


where the last c-number term accounts for the C-function regularized zero-point energy of the 
excited modes: 


TT M 
2L~2' 


^ 2 n- 


TT N-M 
2L 2 


E /r, , TT 3M-N -.x TT iV-3M 

( 2 n+l) = ^^^C(-l) = 


n—1 

The Hilbert space vectors 

|zt) = 


n—0 


... a 


a 


... a 


L 48 


|fci... Um) 


with ji, ■ ■ ■, ji < M and j[, ■ ■ ■ ,j'ii > M form a basis of eigen-states of H with 

H\u) = ^ (2 (fe, T-ifc) + n+ Izz) , 


(4.49) 


(4.50) 


(4.51) 
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where 


= E 2 < + E(24'+1)- 


(4.52) 






The energy density and the energy current in the wires correspond, respectively, to operators 
Ki{t,x) = T[{t,x) +Tl{t,x) and Kl{t,x) = T[{t,x) - Tl{t,x), where 








T:{t,x) = ^^-.Jl{t,xf:- 


A8L^ 


P 

J- 11 


96L2 ’ 


Pd 


961/2 ii ; 


(4.53) 

(4.54) 


are the left-moving and right-moving energy-momentum-tensor components. The constant terms 
are the zero-point energy contributions. Note that the above choice assures by virtue of relation 


N 

'Yd Pii = M that 

i=l 


where 


N 




2=1 


1^ 

Hi{t) = J K'd{t,x) dx 


(4.55) 


(4.56) 


are the observables representing energy in individual wires. 


5 Equilibrium state. 


The equilibrium state at inverse temperature j3 and (electric) potential V is described by the 
density matrix 


PP,v 



Q-m-vQ) 


(5.1) 


where Zjsy is the the partition function. Note that with our conventions, positive V plays the 
role of a positive chemical potential for electrons and of a negative one for holes. Z^y is easily 
calculable with the use of relations (4.51) and (4.39): 




= e 


(N-3M)7: 

48Z 


E- 


-^(k,T-^k)+l3V{p,T-^k) 


' Vein) e' 


■5r^ 


n>0 


M 


2^P' 

n>0 


.(n) e’ 


■5r^ 


N-M 


(5.2) 


with Pein) (poin)) standing for the number of partitions of n into a sum of even (odd) numbers. 
The Poisson resummation formula applied to the fc-sum and the standard relation of the generat- 

OO 

ing function for partitions to the Dedekind function ? 7 (r) = e^ 0(1” allow to rewrite 

n=l 

(5.2) as 


Zpy 


M 


L0V^ Y- „2 
e 47r ^i'i 


E [,(§)]»-“'[,(f(5.3) 


or, with the use of the modular property of ? 7 (r) = ^ ^ ■^ vi~i 

feGZ" ^ 


The equilibrium state uj^ y 
defined by the formula 


expectations of the observable algebra generated by currents 


Jf 


are 


‘^p,vi^) — Tr(p/3,y4) . 


(5.5) 
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The superscript L in wj" y stresses that the state pertains to the junction of wires of length L. In 
forming observables, it is enough to consider only the currents J^{t,x) at fixed t and real x. We 
shall decompose such currents into the contributions from the zero modes and the excited modes: 

+J-{t,x), (5.6) 


see (4.23). In the equilibrium state expectation of products of currents, the contributions from 
the zero modes and from the excited modes factorize. In particular. 




.v(n 


2 L *^04 I 


fc=l 


{2pVL)K 


K 

n 

fc=i 


{2pVL)K 


t fc=l m fcez“ 

g-^(fc,T-ife)+/3y(p,T-ifc) 

feez" 

feez»^ 

fcez'^ 


n « E , 


, (5.7) 


where the second equality results from the Poisson resummation. On the other hand, the expec¬ 
tations of products of are calculated by the Wick rule with 

uj^^y{jfit,x)) =0, (5.8) 

^^,v{Jtit,xi)Jl{t,X 2 )) = fe{xi -X 2 ) + - P,y^)fo(xi -Xa)), (5.9) 

where 


fe{xi - X 2 ) = p{xi - X 2 ; 2L, -i/3) -I- Ce 

fo{xi - X 2 ) = 2p{xi - X 2 ; 4L, -i/3) - p(xi - X 2 ; 2L, -i/3) -f Co (5.10) 


with the constants 


n^O 

a =(i)X -1 - 11 ; + 1 : »»h-^(ii))- 


n^O n^O 

Above, p{z] 0 Ji^ijj 2 ) = p(z;uj 2 , 0 Ji) is the Weierstrass function of period cci and C 02 [47]: 

1 1 


p(z;wi,W 2 ) = ^ + ( 


n^+m^^O 




(z+mcji+ncj2)^ (muji+nuj2) 

— 2 / TT^z—nuJi) 


n^O 


(5.11) 

(5.12) 


(5.13) 


Note the singularity of the 2-point functions (5.9) at the insertion points coinciding modulo 2L. 
For such points, the equal-time commutators of currents have contact terms, see (4.33). For the 
1 -point function of the left current, one obtains 






7n ^—\LV {p,h) — ^(k^Tk) 


^-iLV(p,k)-^{k,Tk) 

fcez" 


(5.14) 


where we have used the relation 

■EcA^(^(p,t-v)) = E E^*. = t • 


2L/3y 


(5.15) 
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From (4.32) or (4.30), it follows that so that 

y (jf(t,a:)) , y (j/(t, x)) = 0. (5.16) 

Hence, in the equilibrium state, the mean charge density is constant in each wire, whereas the 
mean current vanishes. 

The equilibrium state y is invariant under the replacement -O- , the property express¬ 

ing its time-reversal invariance. For the energy-momentum tensor components Tf’^{t,x) defined 
by (4.53) and (4.54), we obtain: 

w^yy{T^{t,x)) = + (1 - 

= jRt,x + e) J-(t,x) + 

= x + e)Jl {t, x) + = uj^y {T[ (t, x)), (5.17) 

which is a consequence of the operator product expansion 

^ Jf (t, x + e) Jf{t, x) = -j^+ x) + ... (5.18) 

^Jl{t,x + t)Jl{t,x) = +Tl{t,x) + ... (5.19) 

holding under the equilibrium expectations away from other insertions points. 


6 Functional integral representation 

6.1 Case with V = 0 

For H = 0, the partition function Zp^Q = Zp and the expectations in the thermal equilibrium 
states ujp g = LOp may be represented by Euclidean functional integrals over a cylindrical open¬ 
string worldsheet, see Fig. 2. 



Figure 2: Open string worldsheet 


For the partition function, 

Z/3= j e-‘^®[slP5 (6.1) 

where the functional integral is over the maps g{t,x) = from ffi. x [0,L] to 17(1)'^ 

periodic in t 

g{t + I3,x) = g{t,x) (6.2) 

with the boundary conditions 


g(t,0)GH, P{g ^d^g){t,Q) = Q , {g ^d^g){t,L) =Q 
and the Euclidean action functional 

N P L 


(6.3) 


5'is[£/] = 4^^ / + (^^‘7i)^)(i:2:)dx = . (6.4) 


0 0 
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To give sense to the functional integrals, one decomposes the multivalued fields ipi(t,x) into the 
linear part which winds in the time direction and the periodic part: 

ip^{t,x) =, (6.5) 

where 

M M 

ni='^K^qjn, = '^ , P da;(p{t, 0) = 0 = da;(p{t, L) (6.6) 

m—1 m—1 


with Qm G Z, iprn G K, and with the multivaluedness reduced to that of defined modulo 2'k. 
The Euclidean action functional decomposes accordingly: 

SeIv] = ^ (fc, Tk) + SsiP^) + Se{{I - P)^) (6.7) 

leading to the factorization of the functional integral 

^ /e"'^®[^‘^'T>(P(p) / P((/- P)<p). (6.8) 

■’ fcGZ" 

The last factor is a standard Gaussian functional integral with the quadratic form corresponding 
to the Laplacian with the periodic boundary conditions in the t direction and the mixed Dirichlet 
one at a; = 0 and the Neumann one at a; = L in the x direction. Such Laplacian is strictly 
positive. Using the zeta-function regularization of such an infinite-dimensional Gaussian integral, 
one obtains: 




V{ii-P)<p)= v{^) 


N-M I -(N-M) 

-idr) 


In the first functional integral on the right hand side of (6.8), we parameterize 


M 


P(p=^ 


(6.9) 


( 6 . 10 ) 


Se[P<p] becomes then a quadratic form in (f/'m) corresponding to the Laplacian with the periodic 
boundary conditions in the t direction and the Neumann ones in the x direction, with constant 
zero modes. The zero-mode integration may be turned to a one over U{1)^ using collective 
coordinates and recalling that fields are determined modulo 2tt. Employing the zeta-function 
regularization for the remaining Gaussian functional integral over the other modes, one obtains 


_ . . M 

e-spP^]v{P<^) = V^[^y \v{^) 


-M 


( 6 . 11 ) 


Upon the substitution of (6.9) and (6.11) to (6.7), the functional integral expression (6.1) for Zp 
reduces to (5.3) with U = 0. 

The expectations of products of equal-time currents in the thermal state are represented 
by the normalized functional integrals: 

P = ^ / Vg , (6.12) 

fc=l fc' = l ^ fc k' 

where on the right hand side 

j^it,x) = jl{t,x) = + idt)(pi{t,x), (6.13) 


are functionals of field ip{t,x) that in terms of decomposition (6.5) take the form 


= -i^ + - idt)^i{t,x) 


277 2 ' 


jKt,x) = -i^ - ^^id^ + idt)‘fiit,x). (6.14) 


The functional integral (6.12) factorizes similarly as in (6.8), with terms —X) 

m 

contributing to the factor with the sum over k and terms with derivatives of (fi entering the 
factors involving the Gaussian integrals calculated by the Wick rule. The latter leads to combina¬ 
tions of products of derivatives of the Green functions of the Laplacians that reduce to expressions 
involving the Weierstrass functions. At the end, one obtains the same formulae as the U —>■ 0 
limit of the ones worked out before for the expectations of products of the left-moving currents 
resulting from applying the rule (4.30) to the right-moving currents. 
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6.2 General case 


An imaginary potential V may be included in the functional integral approach by imposing the 
twisted-periodic boundary conditions in the time direction on the C/(l)^-valued fields g = [gi) = 

(eivi). 

gi{t +I3,x) = . (6.15) 

The latter may be implemented in the functional integral by decomposing 

LPi{t,x) = {iV +‘^ni)t + (pi{t,x) (6.16) 

with ipi periodic in the time direction, keeping the same boundary conditions in the x direction 
that take again the form (6.6). For real V, the above decomposition implies a complex shift of 
the functional integration contour over fields g. Performing the functional integration the same 
way as before, one obtains the representation 

Zg.v = J Vg, 

k—1 k' — l 

where the currents are still given by Eq. (6.13) 
on V in the way described by decomposition 

7 Closed-string picture 

7.1 Classical description 

A symmetric role of time and space in the functional integration leads, upon reversing those roles, 
to a description of the equilibrium expectations in the closed-string picture, see Fig. 3. 


(6.17) 


/ (6.18) 

k k' 

and the contour of functional integration depends 
(6.16). 



Figure 3: Closed string worldsheet 


In the latter, a collection of N closed strings of length /3, is described by fields gi{t, x) = 
defined for real t and x and twisted-periodic in the x direction: 

g^{t,x +P) = g^{t,x)e~^^ , (7.1) 

where V is taken imaginary, compare to (6.15). On the classical level and for Minkowski time, 
such fields are governed by the action functional 


S[g] 


-T 

477 / ^ 


P 

dt J rmdt(pi)^ - dx. 
0 


The twist in the periodicity condition may be absorbed by setting 


(Pi{t,x) = (fii{t,x) +\Vx, 


(7.2) 


(7.3) 


where 

(p^{t,x + p) = (f)^{t,x)+ 27:1111, rmeZ. (7.4) 
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The classical solutions have the form 


^i(t, x) = ipi{t + x) + (pi (t-x), 


(7.5) 


where 


r/j. I \ ^i.r , y/^TT ^-.T (i I \ I i \ '' 1 

’ (t ± a:) = (Po’ + ± 2:) + ^ ^ - 


az: e 


n^O 


for Oirf = <^o + (^0 = e (R/27rZ)^, + ag) = po ^ and 

1 / / 7- \ 

^(tto - tto) = m 


(7.6) 


(7.7) 


where m is the vector of N winding numbers rrii G Z. The symplectic form on the space of 
classical solutions is 

n ^Spo- ASipo - I Z J ■ ^"(-n) ■ ("^-S) 

n/0 


7.2 Quantization 

The Poisson brackets obtained from ft lead to the following canonical commutators: 

[T’OijPOj] = i^i ^ij J [o^ni > '^n'j ] ~ ^ij'>T‘^n+n',0 t = Tj Sij U Sn+n',0 ■ ("^'9) 

For fixed winding numbers, the zero modes will be represented in the Hilbert space L^(U{lZ) 
with an orthonormal-basis vectors 

N 

\k) = \Z ... k^) = exp (I'y^ k^poi) for k'' £ Z (7-10) 

i=l 

such that 

N 

Zpq, |fc) = -i ^ exp (i ^ fc*Vo*') = fc* |fc) • (7.11) 

The Hilbert space of states for the zero modes is a direct sum of an infinite number of copies of 
LZu{ir), one for each winding vector, 

Ho = © Z{U{\Z), (7.12) 

with an orthonormal-basis vectors |fc,Tn). The non-zero modes are represented in the tensor 
product of two standard Fock spaces generated by applying products of the with negative 
n to the normalized vectors |0)^’’', annihilated by with positive n. The scalar products are 
defined by demanding that (a^’’')l = The Hilbert space of the full theory is 

n=Ho^F^^F^ (7.13) 

and we identify |fe, m) = \k, m) © | 0 )^ © | 0 )’'. 


7.3 Current, energy and (magnetic) charge 

As before, we define the left and right current for the closed string 


jritF) = 


Ill 

2tt 2 


{dt ± dZpiit, = 


7rin(t±3;) ir'^Y 

i H— 

477 


The classical Hamiltonian of the system is 


13 

^ i + {d.Pif + 2iVd,pi - H") dx . 

^ 0 


(7.14) 


(7.15) 
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Once quantized, its t^-independent part becomes the standard Hamiltonian of N closed strings 




where 


j-e,r _ rf 
^Oi - 2 ■ 


l,r 

^ni ^—ni 


(7.16) 

(7.17) 


and the term comes from the zero-point energy. In the action on I-Lq vectors, 


Hr^ \k. m) = 


^ + 0 ^w 2 ) - \k,m) . 


(7.18) 


The action of excited mode operators for positive n raises the eigenvalue of Hcs by 

The part of the Hamiltonian linear in V is equal to iH Q^, where 

QTs = J2 f J[{x,t))dx= ( 7 - 19 ) 

i -fo i 

is the total magnetic charge of the closed (untwisted) strings. It acts only on TLq: 

Q^^\k,m) = Y^r‘^m^\k,m). (7.20) 

I 


Finally, the part of the Hamiltonian quadratic in V is an additive constant, so that the full 
quantum Hamiltonian of the closed-string system becomes 

i7 = i7es + il"QS (7-21) 


7.4 Boundary states 

In the closed-string pictures, the boundary conditions in the space direction, which in that picture 
becomes the time direction, are represented by the boundary states in the (completion of) the 
closed-string space of states [38, 21]. The boundary state that corresponds to Neumann boundary 
condition for all field component is 


\W)) = A^ 


E 

mez'^ 


- E 


|0,m) 


(7.22) 


where Hjv is a suitable normalization constant. This boundary state satisfies the relation 


9t¥>,(0,x)||V)) = 0 (7.23) 

whose excited-mode part implies that 

(aL+y_„)J||Ar))=0 (7.24) 

determining the form of the Ishibashi-type dependence of HA/”)) on those modes. For the mixed 
Dirichlet-Neumann boundary condition describing the junction of wires, the boundary state has 
a more complicated form 


||H)) = Adn (27r)“-^Vd^ 


< E E 

fcez'^ 

mi=J2iK’P's„ 


-EE 




+ E E 


j = M + l n = l 




\k,m), (7.25) 


where Sm run through integers and = A, • a^'^, see (4.9). One has 

PdMO, 11^)) =0, (1 - P) <p(0, x) I \B)) = 0, 


(7.26) 
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where field tp may be equivalently replaced by (p. The excited-mode part of these conditions 
implies that 


(«l,+fi(-„b)ll^)) = 0 for j<M, (d^-d[_„)^.)ll^)) = 0 for j>M, (7.27) 

fixing the form of the Ishibashi building-blocks of | \B)). The zero-mode part of the first of relations 
(7.26) assures that ^ = 0, whereas the zero-mode part of the second relations implies that 

i 

(1 — P)m = 0 which is solved by ^ for integer Sm- The sum 

m 

(27r)^-^ Vdetr ^ |fe) (7.28) 

fcGZ" 

fe‘Kr=o 

represents the delta-function supported by the brane B = k{U{\)^) C 17(1)^ defined by the 
integral 

^b{po) = y it ^ivoi - Vdet T J] dj^rn (7.29) 

i m m 

over U{\)^ of the 27r—periodic iV-dimensional 5-function. Indeed, 

fcGZ'^ m 

= ^ (7.30) 


which reproduces (7.28). 


7.5 Partition function 


In the closed-string picture, the partition function Zpy is represented by the matrix element of 
the Euclidean evolution operator between the boundary states. A direct calculation gives: 


Z 0 ,v = ((A/’l I e 

= An\/ detT |??(^) 


/ 2iL \ 

N- 2 M 

1 




m 

-)J 




E 

i=Yl 


3-^ Y,{rimiY-iLV 


(7.31) 


This coincides with expression (5.4) upon relabeling s = k and recalling the definition (4.40) of 
vector p, provided that 

AnAb{ 2Tr)^-^ = 2-^ . (7.32) 

The latter identity is assured if we take An = A^ and Ab = A^for A = -^ and 
B = 


7.6 Expectations 

The expectation values of products of currents in equilibrium state y take in the open string 
picture the form of the matrix elements between the boundary states of the time ordered products 
of Euclidean versions of currents 
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where the Euclidean time ordering puts the operators at bigger Xk or yk' to the left. The 
powers of —i and i represent the derivatives of the Euclidean conformal change of variables, 
X + it 1 -^ t — ix and y — it t + iy, respectively, that reverses the roles of time and space. 

The proof of (7.33) in done in few steps. First, consider only the left currents. As in the 
initial picture, we distinguish the constant part from the excited terms, 


Jl{t,x) 




(7.34) 


see (7.14). This decomposition factorizes in the expectation values. For the constant terms, we 
get by direct calculation: 
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(7.35) 


fc=i 




which almost reproduces the zero mode part expectation value (5.7) of the initial calculation but 
with one extra term, where {... (fcp. Ip )...} runs through all possible pairing of {1,..., K}, as in 
the Wick theorem. The presence of this term can be seen by induction on K. On the other hand, 
the expectations of products of are calculated by the Wick rule with 

- iZ^^y {{Af\\e-^^jf{-ix, 0) m = 0, (7.36) 

i-i)^Zp^y m I (-ixi, 0)Jl {-ix2 ,0) 11^)) (7.37) 

= “I (^) (^*1*2 {p{xi - X2] 2L, -i/3) + C') 

+ {6iy^ - Piy^){2p{xi - X2;4:L, -i/3) - p{xi - X2\2L, -i/3) + C')^, (7.38) 


where we get expressions with the Weierstrass function similar to (5.9) but with different constants 

p^O 

Co = {^) (|-2^[sinh(7ri|^)]-" + ^[sinh(^^)]-^). (7.39) 

p/O p^O 

The theory of Weierstrass function of periods uji and 1 x 2 [47] provides the identity 

S(1 + E ("i^)) - /rd + E (^)) = ±i. (r40) 

n^O n^O 

where the sign on the right hand side is that of the imaginary part of a;i/w 2 - This leads to the 
relations 


C'; = C'o, C'=Ce-^. (7.41) 

The contribution from the last term will cancel exactly the last contribution appearing in (7.35) 
establishing identity (7.33) for any product of left currents. Finally, the closed-string expectation 
value of a general product of left and right current will be a combination of factors corresponding 
to the decomposition (7.14). By direct calculation, the constant part of the right currents can 
be expressed via the S matrix in terms of the one of the left currents with the use of (4.26) and 
the fact that S also preserves vectors kP. In the computation of the excited part, the S matrix 
appears naturally upon noticing that in the proper basis defined in (4.9), it becomes 

^ = diag(d_^,-!,...,-!) (7.42) 
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which is precisely how the excited inodes and are related when they act on ||yB)), see 
(7.27). Finally, under the closed-string expectation every right current is related to the left one 
by the S matrix, exactly as in the initial picture (4.30). This proves identity (7.33) in the general 
case. 


8 Thermodynamic limit 

In the thermodynamic limit L —>■ c» the wires become infinitely long. The partition function 
Zfj^v diverges in that situation but the free energy per unit length has a limit: 
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( 8 . 1 ) 


as easily follows from its form (5.4). The equilibrium state expectation values of the products of 
currents also possess the L —>■ oo limit. In particular, it follows from (5.7) that 
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for real V and relations (5.8) and (5.9) imply that 
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( 8 . 2 ) 

(8.3) 

(8.4) 


as both fe(x) and fo{x) tend to sinh ^ when L ^ oo. The latter property follows 

from (5.13) and the identity (7.40). Eqs. (8.2), (5.8), (8.4) and the Wick rule, as well as the 
relation (4.30), determine the L —>■ oo limit of the states ojj^y. Unlike for finite L, that 
limit is not represented by a trace with a density matrix (for L = oo, the Hamiltonian has a 
continuous spectrum and the operator is not traceclass). In particular, one obtains: 


w/3,y(jf(t,a;)) = 
i^0yiJt(t,xi)Jl(t,X2)) 
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The operator product expressions (5.17) and the limit (8.4) (that is uniform in small |a;i — X 2 \) 
imply that 


w/3,y(7’/(t,a;)) 
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(8.7) 


In particular, the mean energy density in the equilibrium state is constant in each semi-infinite 
wire (but differs from one wire to another) and the mean energy current vanishes. 

The L = 00 state is easy to represent in the closed-string picture: by examining the right 
hand side of (7.33), one infers that the boundary state \\Af)) of (7.22) is projected when L —>■ 00 
to the closed-string vacuum |0, 0) so that 
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for Xk,yk' > 0 , 

One of the crucial observations that follows from (8.2) and (8.4) is that, when restricted to 
the products of left-moving currents Jf{0,x) with a; > 0, the limiting L = 00 equilibrium 


23 



expectations do not depend on the choice of the brane B describing the contact of wires. In 
particular, such expectations are the same as for the space-filling brane Bq with Sw = 5^ 
corresponding to the disconnected wires for which Jl{t,x) = Jf{t,—x) and ||,Bo)) = IIA/”)). The 
physical reason for this behavior of the expectations of left-moving currents is that the latter 
did not have contact with the junction up to time zero. The above observation is essential 
for the construction of nonequilibrium stationary state where the individual wires are kept at 
different temperatures and at different potentials. For the disconnected wires, one has the obvious 
factorization: 

K N 

jHt,Xk)y ( 8 . 9 ) 

k—1 i—1 k 

Hence the same formula holds in the L —>■ c» limit of the equilibrium state for any brane B = 
For the disconnected wires, one also has the relation: 

K' N N 

n= nnn(s.io) 

fc' = l i=k' k' *=1 k' 

iy=i 

It is easy to check using (4.30) and (4.26) that the latter factorization holds in the limit L —>■ oo 
also for other branes B = 

9 Nonequilibrium stationary state 

Following [32, 11], see also [41], we shall consider a nonequilibrium stationary state (NESS) Wneq 
describing the situation when different semi-infinite wires are kept at different temperatures and 
different potentials. State Wneq may be obtained by the following limiting procedure. For each 
disconnected semi-infinite wire, one considers the algebra Ai generated by products of currents 
j/’’'(0,x) for X > 0, together with a state given by the restriction to Ai of the L = oo 

N 

equilibrium state 0 J| 3 ^y^ for the space-filling brane Bq. The product state Win = (S' y. on 

algebra A = 0i Ai describes the disconnected wires with each prepared in its own equilibrium 
state y.. As in Sec. 4, algebra A may be identified with the one generated by currents 

J^’^{0,x) with X > 0 and 1 < i < N for the connected wires. Let I4t for t > 0 describe the 
forward in time Heisenberg-picture evolution of the currents J^’^{0,x) with x > 0 in the presence 
of brane B: 
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Jn 0 ,x-f)) 


Jn0,x-f) 

Y.S^^'4{GA 

i' 


for t < X 
x) for X <t 


(9.1) 

(9.2) 


Then for A ^ 
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In order to prove the above relation consider the backward in time Heisenberg evolution for 
decoupled wires, i.e. in the presence of brane Bq: 


W°,j/(0,x) = J/(-Lx) = J/(0,x- f) = 


jno,t-x) 


w°,jno,x) = jr(-t,x) = jr(o,x+t) = j:(o,t+x) 


for X, t > 0. Such a decoupled evolution preserves the product state 
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for X <t, 

Win so that 
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is the scattering operator in the action on algebra A. The explicit form of the latter in the action 
on the chiral currents follows from equations (9.1), (9.2), (9.4) and (9.5): 

sj!(0,x) = Jf(0,x), 

SJ[(0,x) = 

i' 

for a; > 0. Note that the nonequilibrium state Wneq is preserved by the Heisenberg evolution Lit 
so that its stationarity follows. Hence the explicit formula: 
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(9.9) 
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where on the left hand side the currents correspond to connected wires and on the right hand side 
to disconnected ones and the values of t, Xk and Xk' may be taken arbitrary (with noncoincident 
xi,..., Xk, —yi, ■ ■ ■, —yK' to avoid singularities). In particular. 




K 


k=l 


n aa^>a;fc) ) = Win! 


K 


k=l 


(9.11) 


so that the difference between Wneq and Win, due to the junction between wires, arises only in 
the presence of right-moving currents. The left-moving currents do not feel the influence of the 
junction. It should be stressed that the dynamics considered above both in the presence of 
the junction and for decoupled wires is generated by the Hamiltonians that do not include the 
electric potentials in the bulk of the wires. Those play the role only in the preparation of the 
initial product state and may be applied far away from the junction. That the ballistic evolution 
of chiral currents persists for long times in the bulk of the wires in such a nonequilibrium situation 
should be assured by the integrable nature of the Luttinger liquids, see the discussion at the end 
of [13]. 

Specifying Eq. (9.10) to the 1-point expectations, one obtains: 
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so that the mean charge density and mean current in the wires are 
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(9.13) 
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respectively. They are constant in each wire and the mean current does not vanish, in general, 
at difference with the equilibrium state. The electric conductance tensor of the junction (in the 
units e^/h) is 


el^ — 9 , , 

— dV / 
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(9.14) 


This agrees with the calculation of [39, 40] based on the combination of the Green-Kubo formula 
with the conformal field theory representation of the equilibrium state. Note that the conductance 
vanishes for the decoupled wires. The nonequilibrium current 2-point functions are given by 
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Note that the nonequilibrium states Wneq with coupled wires break the time reversal symmetry 

For the expectation value of the energy-momentum components, we obtain from the operator 
product expansions (5.18) and (5.19) 
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SO that the mean energy density and current are, respectively, 
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This results in the thermal conductance 
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10 Full counting statistics 


10.1 Charge transport 

Measuring transport of charges through the junction of quantum wires requires specifying mea¬ 
surement protocol that may be not easy to implement. Refs. [29, 30] proposed an indirect mea¬ 
surement of charge transferred through a quantum resistor and obtained a closed formula for 
statistics of the results. The same charge transfer statistics could be obtained by considering a 
direct two-times measurement of the total charge accumulated in the system, provided the latter 
is finite. Following [11], we shall employ the second measurement protocol that is conceptually 
simpler although unpractical for large systems, keeping in mind that the charge transfer statistics 
obtained this way may be also accessed by a more practical indirect measurement protocol. 

Consider first the system of disconnected wires of length T, each with Hamiltonian and 
charge operator Q^. Prepare the system in the product state Wq = 0 given by the 

density matrix po = <8 Pa V) where v = ^ zero, measure the 

total charge in each wire. Then connect the wires instantaneously and let the system evolve. 
At time t, disconnect the wires and measure the total charge Qiit) in each wire. By spectral 
decomposition, 

Q? = E = E w • (10-1) 

The probability that the first measurement gives the values of charges {q^) = is equal to 
tr((g)A Pp. Vi^iqo)- After the first measurement, the density matrix is reduced to 
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( 10 . 2 ) 


The probability that the second measurement gives the values of charges (pi) = q, is then equal 
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to tr(/?o+ n Altogether, the joint probability of the results {q^^q) is 
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where to obtain the second equality, we used the fact that commute with p'‘^, y.. The 

probability that the charges change by Aqi = Qi — is 
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(10.4) 


The latter probabilities may be encoded in their characteristic function called the generating 
function of full counting statistics (PCS) for the electric charge transfers: 
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For connected wires, the change of the wire charges in time is 

t t L 


(10.5) 


AQi{t) = Qi{t) - Qi{0) = J ■^Q^{s)ds = J ds J ds{J^{s,x) + Jl{s,x)) dx 

0 0 0 
t L t 

= y*ds J {dj,jf{s,x)-d:rJl{s,x))dx = -J (j/(s, 0) - J[(5,0)) ds 

0 0 0 
t 

= - j (j/(0,s)-^5.,.4(0,s)) ds. 


( 10 . 6 ) 


After disconnecting the wires at time t < L, the latter observables become the ones for uncon¬ 
nected wires given by the right hand side of (10.6). The crucial fact is that they are extensive in 
time but not in the wire length, unlike the total charges. Note the commutation: 
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[AQ,(t),Q0] =- j ds j [(J/(0,s)-^5,,»4,(0,s)),4(0,x)] dx 
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= —{Sii' — Sii') J ds J (5'(s — X -f 2nL) dx = 0 . 


(10.7) 


Since the observable Qi{t) become equal to Q® -I- AQiit) after the disconnection of wires, the 
PCS generating function (10.5) may be rewritten due to (10.7) in the simpler form 
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Due to the translation invariance of the state 7 
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In the limit L —>■ c», the initial states Wg tend to the product state Win for semi-infinite wires 
considered in Sec. 9 so that 
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with the last equality following from relations (9.11) and (10.6). 

We would like to study the large deviation form of the FCS generating function by calculating 
the rate function 

( 10 . 12 ) 


7(i/) = lim y \n^'Ft{u) . 


Refs. [11, 19] exposed a strategy for the calculation of such rate functions for semi-infinite wires 
with a purely transmitting junction from its derivatives. Applying it to our case, we note that 
such derivatives have the form 
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with Vi as above. It was argued in [11, 19], following the approach set up in [10], that 
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because for large t the exponential factor becomes close to e * providing effectively the 
imaginary additions to potentials Vi. Since the one point function on the right hand side of 
(10.14) is independent of s', this line of thought gives by the analytic continuation of (8.5) the 
identity 
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which, together with (10.9) and the relation /(O) = 0, implies that 
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The existence of the limit (10.12) means that at long times the PDF of charge transfers takes 


the large-deviations form 

Pi(Aq) ~ 

(10.17) 

where the rate function 
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®'/(p) = max(EPi^» “ ®y(-ii^)) 
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(10.18) 


is the Legendre transform of *^(— 11 /). For given by (10.16), ®*/(p) is a quadratic poly¬ 

nomial on the subspace where it is finite. In other words, the large deviations of charge transfers 
per unit time have the Gaussian distribution with mean 
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equal to the mean current in the nonequilibrium state, see (9.13), and covariance 
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(10.19) 


( 10 . 20 ) 


Note that the first of equalities (10.16) implies that the large-deviations rate function for 
FCS of charge transfers is proportional to the difference of equilibrium free energies for different 
potentials: 

N 

"'/M = (10.21) 

i=l 

where fi{P, V) = —is the equilibrium free energy per unit length in a single decoupled 
semi-infinite wire with Neumann boundary conditions, see (8.1). Relation (10.21) implies in turn 
that 

= hm i (10.22) 

t=2L—>-00 t 

if we define ®^F’/^'(i>') for t > T by the right hand side of (10.8). Indeed, in that case the 2"'^ 
equality in (10.10) implies that 
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where the partition functions on the right hand side pertain to the disconnected wires of length 
L. A priori, it is not clear that the same result for ^'/(ly) arises in the physically different limit 
that takes the thermodynamic limit L —>■ oo before sending t —>■ oo. The calculation of [11, 19] 
amounts to the claim that both limits are equal. 


10.2 


Exact result for 




The exactly soluble nature of the model considered here allows to examine closer the distribution 
of charge transfers for finite L and t and to see in more details how its large-deviations form 
arises. A direct calculation performed in Appendix A gives the result 
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where the subscript “reg” refers to a necessary ultraviolet regularization, that replaces the diver¬ 
gent constant Coo = Z) ^ by 

n>0 


AL 

CAi = E^=MAC) + C+ O(^) (10.25) 

n—1 

with the ultraviolet cutoff A, see Appendix A. Variables are as before, see (10.9), and 

6»i(r; z) = J2 e’^‘^(”+5)'+2’^*("+5)C+5) (10.26) 


is one of the Jacobi theta-functions. The first exponential factor on the right hand side of (10.24) 
is the characteristic function of a centered Gaussian distribution of charge transfers Aq with the 
covariance 
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The t-dependent expression under the logarithm is positive for 0 < t < 2L. Note that the 
ultraviolet divergent contribution to the covariance is independent of t. It describes the charge 
transfers that arise at the moments of the connection of wires at time 0 or their disconnection 
at time t but do not contribute to the average charge transfers realized during the long period 
of time when the wires are connected. The second factor on the right hand side of (10.5) is a 
characteristic function of the discrete distribution 


exp[j2 U 
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E exp[ E 
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of charge transfers Aq. The two types of charge transfers are realized independently and both 
correspond to the vanishing total charge transfer As we shall see below, they both 

contribute to the large deviations result (10.16). 

In order to study the behavior of the charge-transfer distribution for large L and large t, we 
shall rewrite (10.24) applying the Poisson resummation formula to the fc-sums and the modular 

transformation 9 i{t; z) = i(—ir)”^/^ e“^^ 0i(—f-) to the Jacobi theta function. The resulting 
expression is 
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Together with relation (10.25), it allows to extract the large L behavior 
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where c > 0 is some Pi- and rj-dependent constant. The first exponential factor describes the 
leading behavior of the contribution in the 1®* line of (10.29) and the second one that in the 2"'^ 
line. We infer that 
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and 
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reproducing the large deviations result (10.16) up to the ultraviolet regularization. Note that 
if follows from relation (10.30) that the same result is obtained for the limit of ln®^i^/'(i/)reg 
obtained by sending simultaneously A, L and t to inhnity in such a way that the ratios and 
tend to zero. This specifies more precisely the region where the distribution of charge transfers 
takes the Gaussian large deviation form (10.17) described previously. The above analysis does 
not cover, however, the case (10.22) with t = 2L —>■ oo which, although giving the same limit, is 
somewhat special. In particular, no ultraviolet regularization of is required for ®^T'^(t'). 


10.3 Heat transport 

The protocol for the measurement of the thermal transfers is the same. It consists of preparing 
the system of wires of length L in the initial product state ojq = performing 

the measurements of the energies Hi (0) = and Hi (t) in the disconnected wires at two times 
in between which the wires were connected. Denoting the results, respectively, (e°) = e° and 
(ci) = e, we encode the probability of the change of energies Ae^ = 6^ — 6° in the characteristic 
function 

Ae 

the generating function of FCS for heat transfers. The change of energy of the wires connected 
between times 0 and t is 
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AHi(t) = Hi(t)-Hi{0)= J ^H^{s)ds = J ds J ds(Tf{s, x) + T[{s, x)) dx 

0 0 0 
t L t 

= JdsJ {d,T^{s, x) - d,T[is, x)) dx = - y' (1)^(0, s) - T[{s, 0)) ds , (10.34) 

0 0 0 

compare to (10.6). Moreover 

T[{s, 0) = T[{0, -s) = lim (I J[(0, -s + e) J[(0, -s) + 

= JTo (f ^ SieSu"4{0, s - e) 4,(0, s) + 

i' i" 

= ^Y.4TT’frlTl,{d,s) + 4, Y. Su'44s)Su"44,s) (10.35) 

i' i' 


so that 


= /(5]4*'-4»')")41^/'(0,s)- 2^^ ,Si,,4(0,s),S,,„4,(0,s))ds. (10.36) 

* do j, i'/i" 

Interpreting the latter operators as observables for disconnected wires, we have the commutation 
relation 


[Hl/^H,{t)\ 


t 

^ J{4j - : (4(0,s)) 2 : -2 ^ ^,,,^,,4 (0, s) 9,4(0, s)) ds 

{44t)? ■■ - ■■ (44^4 ■■) 

t 

- ^ E / 4 (0, s) 9,4(0, s) ds 

0 

(10.37) 
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as a consequence of the identity 


[H°,jf(t,x)] = -iS^jdtJ-(t,x) = -iSijda;J-{t,x). 


(10.38) 


Note that [H^, AHi{t)] ^ 0 and the generating function (10.33) of PCS for heat transfers 


th^L 


(A) = 




). (10.39) 


This difference occurs even for t = 2L when the first term on the right hand side of (10.37) 
vanishes, but not the second one. 

We shall calculate explicitly *^F/"(A) for t = 2L which is easier than for general t. In this 
case, 


77° + AH,{2L) = Oii'Oil" ^O^QiCtQi' ^ ^ ^ Ct—2n,i^2n,'. 

i'%" n>0 


(10.40) 


in terms of the modes, where O is the orthogonal matrix related to matrix S by (4.27). One easily 
checks that the above observables commute so that they may indeed be measured simultaneously 
in the disconnected wires. Let 

^ ^ A*(77° + A77,(2L)) = ^ ^(OAO)*,. + 2 ^ a.2n,id2„y) , (10.41) 




n>0 


where A stands for the diagonal N x N matrix with entries Ai. The contributions of the zero 
modes and of the excited modes to the expectation 




-J2 0iiH°-VQ°) , 


iyt 


Wn e 


factorize. The first one has the form 


tr« (^1 


-T,Pii.H°-VQ°) 


(10.42) 


^ g-f (r“^fc,,3r“^fc) + (/3V,fe)- , Cr“^fe) 

fcez" 




= det {l + ip-^C) 


0— ^ 5 /3r“^fc) + {^V , k) 

^ {rk , (P+iC)-^rk) - iL (rjSV , (/3+iC)-Vfc) + i(r/3V , (,3+iC)~ V/3V) 

fcez" 

^—ttL (rk , l3~^rk) — iL {rl3V , I3~^rk) +-^{rlSV , 0~^rl3V) 

fcez" 


(10.43) 


where r and /3 stand for the diagonal matrices with entries (r^) and (Pi), respectively. 


C = X- 0X0 (10.44) 

is a symmetric N x N matrix, and V denotes the vector with components Vi. The right hand 
side of (10.43) was obtained by the Poisson resummation. As for the contribution of the excited 
modes, its calculation is given in Appendix B and results in 

det ( 7 + ( 7 -e"T">' 0 eT"^ 0 )(eT /3 _ (10.45) 

n>0 

with a convergent infinite product. Gathering expressions (10.43) and (10.45), we obtain: 

= det (I + il3-^Cy^^'^ 

{rk , (/3+iC)-Vfe) - iL {rfiV , (/S+iC)" Vfc) + ^(r/3V , (/3+iC)-T/3V) 

fcez" 

^ g— ttL (rfc ,/3“irfc) — iL (r/5V ,/3“irfc) + ^(r/3V ,/5“ir/5V) 

fcez" 

X det ( 7 +(7-e-T"^OeT"^0)(eT/3 _ (10.46) 

n>0 
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In the limit t = 2L —>■ c», 


1 


In 


L—¥00 


^ ({rPV , {P + iCy^rfiV) - {rV , prV]) 

oo 

- y In det (/+(/- O 


0 

_ th 


^ ‘7(A) 



dx 


(10.47) 


for sufficiently small |Ai| so the zero-mode contribution with k = 0 dominates. 

If we calculated the right hand side of (10.39) for t = 2L instead of *^F 2 ^(A), the only 
change would be the replacement of matrix e“T”^0e“^"^0 by e“T"(^“OAO) the last 
line of (10.46). In general, however, matrices A and 0X0 do not commute if O has nondiagonal 
elements. Such a modification would also kill the symmetry (B.28) of the contribution (10.45) 
showed in Appendix B. The difference of resulting expressions would persist also in the L ^ oo 
limit of (10.47). 

An explicit calculation of *^F/"(A) for t ^ 2L is also possible along the lines of Appendix 
B, using the expansion of + AHi{t) in terms of the modes. We expect that the same large- 
deviations rate function (10.47) for energy transfers would result if we sent L to infinity before t, 
as suggested by the analysis of [11], but proving that basing on the exact formula for *^F/"(A) 
requires technical work that we postponed to the future. 


10.4 FCS for charge and heat and fluctuation relations 

The characteristic function of joint measurements of charge and heat transfers is defined as 


F; 


^(«^,A) = 


UiAqi-\-XiAei 


Ft{Aq,Ae) 


Aq.Ae 


= UJn{e 




(10.48) 


For t = 2L, it can be easily computed since there is only a change in the contribution of the zero 
modes with respect to the calculation of Subsec. (10.4). Indeed, 


A) = uJo[e 


-iE iY:\i{H°+AHi(2L)) 

" e * 


(10.49) 


so that the only effect is the change of F to V — i/3 ^0 in the numerators of (10.43). We infer 
that 

= det(/ + z/3-‘C)-'/" 

{rk , {0+iC)~^rk) — iL {r{0V—ii>) , (/S+iC)”^rfc) + ■^{r{0V—ii>) , (/3+iC)“^r(/3V —ii>)) 

^ - 7 - 

g— ttL {rk , 0~^rk) — iL {r0V , 0~^rk) + -^{r^V , 0~^r0'V) 




X ]^det(/-b(/-e-T"^eT”‘^^0)(eT'^-/) 


-1 


(10.50) 


n>0 


with 


— \nFLii^,X) 

ZL/ L—^oo 


^ ({t{I 3V - iD ), (/3 + iC)-L{pV - iz>)) - {rV , /JrF)) 

OO 

- J In det (^I + (I - dx 

0 

= /(«^,A). 


(10.51) 


The large-deviations rate function function (10.51) of FCS for charge and heat transfers satisfies 
the fluctuation relation [2, 12] 


f{u,X) = fi-t^-ipV,-X + if3) 


(10.52) 
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that reflects the time-reversal invariance of the dynamics. The generating function (10.50) does 
not possess, however, the corresponding symmetry which arises only in the t = 2L —>■ c» limit. 
Relation (10.52) is a consequence of the following matrix transformation properties under the 
change {u, A) i—^ {—u — i/3V, —A + i/3): 

i/3 + iC =/3 + iX-iOXO ^ /3 + i{-X + i/3) - iO{X + 1/3)0 = 0{(3 + iC)0, (10.53) 

r{/3V - iO) = r/3V - i(/ - 0)ru ^ r/3V - i{I - 0)r{-u - i/3V) = 0{r{/3V - iu)) (10.54) 

and of the symmetry (B.28) showed in Appendix B. That the same symmetry fails to hold for 
the generating function (10.50) follows from the fact that under the change (i/. A) i—>■ {—v — 
i(3V, —A -f i/3) the sum over vectors k G in the numerator of the middle line of (10.50) is 
transformed into the one over vectors r~^Ork that, in general, do not belong to . 


11 Comparison to Levitov-Lesovik formulae 

In [29], L. S. Levitov and G. B. Lesovik obtained a closed formula for the FCS of charge transfers 
between N free fermionic systems, as those of Sec. 2.2. Such systems are assumed to be initially 
in different equilibrium states and to interact subsequently during a period of time t. Their 
interaction is described by an x unitary mode-dependent matrix St(p) accounting for the 
scattering between the fermions of different systems, see also [30]. The Levitov-Lesovik formula 
for the generating function of charge FCS has the form of a product over the free fermionic modes 
of determinants: 


$,(iv) = ndet(/-/(p)+/(p)e-‘^(^’)‘'§*(p)te‘^(^’)‘'St(p)), (11.1) 

p 

where s{p) is the sign function representing the charge of modes, v is the diagonal N x N matrix 
of coefficients Vi and f{p) that of Fermi functions fi{p) = (e^de(p)-s(p)Fi) _|_ 1)“!, with e{p) 
representing the energy of modes. Upon taking the scattering matrix time and mode independent, 
§t(p) = S, and the linear dispersion relation e{p) = ^\p\ as in Sec. 2.2, and upon aligning the 
time and the size of the system by setting t = 2L, the above generating function leads in the 
rate function 


(j){u) = lim i In'hRiz) = f \n ( det (l — f'^ (e) + f'^ (e) e S^e*^ S)) de 

t=2L —>00 t J ^ ' ' 

0 

00 

+ ^ y'ln(det(J-/-(e) + /-(e)e‘^Ste-‘'^§))de, (11.2) 

0 

where /^(e) are the diagonal matrices with entries f^{e) = (e^X'^TVi) _|_ l)“i. Note that this 
is a different expression than the rate function ^'/(iz) of (10.16) obtained in Sec. 10.1 which is 
quadratic in v and V. For closer comparison, let us extract from (11.2) its leading quadratic 
contribution describing the central-limit Gaussian distribution of charge transfers. In Appendix 
C, we show that 


lim 0^ (f){9 ^v) 
6—^00 




= - ^ E - E +E ^ E 




+ESrEis-ivisxoiv-ESEiS' 




+ ^ E A') E . 


(11.3) 


where g{/3i, Pp) is given by the integral formula (C.IO). The first two lines on the right hand side 
reproduce the rate function (10.16) for the compactihcation radii squared rf = 2 that correspond 
to free fermions if we set Su' = p. The last two lines represent terms not present in the rate 
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function (10.16). Of course, in spite of similarities, the coupling between the free fermions realized 
by the junction of wires with matrix S describing the scattering of the cnrrents at the junction 
is different than that assumed in the Levitov-Lesovik approach, so there is no a priori reason 
for the two systems to lead to the same charge transport statistics. Note also that for arbitrary 
unitary matrix (Sji') the matrix (jSii'p) is not necessarily orthogonal. 

In the particular case when all temperatures are equal Pi = P for i = 1,..., N, the last line 
of (11.3) reduces to 

i/' 3-j' i 3,j' 

^ 3:j' 

if we use relation (C.12) and the unitarity of matrix S, and expression (11.3) reduces to 

i i' 

On the other hand, the rate function (10.16) becomes in this case equal to 

®y(«^) = E(^^i + P~^Vi){v, - y] (11.6) 

i i' 

upon using the orthogonality of matrix S = O. It follows that for equal temperatures, 

if we identify |§y P = Sij, assuming that the latter identification leads to a matrix S with the 
desired properties. In that case, the fluctuations of charge transfers indnced by different electric 
potentials at the same ambient temperature agree in the two setups on the level of the Gaussian 
central limit contributions. One should remark, however, that the scaling limit (11.3) removes 
from the Levitov-Lesovik rate function (11.2) the term linear in V and quadratic in u that is 
responsible for the zero-temperature shot noise given by the Khlus-Lesovik-Biittiker formula [14]. 

There is another relation of the PCS statistics that we have obtained for the junction of 
wires and the Levitov-Lesovik type formulae, this time for the energy transfers. Indeed, the 
contribution (10.45) of the excited modes to the generating function of energy FCS 

coincides with the version of the Levitov-Lesovik formula for N free bosons with the dispersion 
relation e„ = ^ and the interaction described by the scattering matrix § = O. The bosonic 
version of the Levitov-Lesovik formula was obtained in [28]. Its proof in that reference provides 
a more direct way to calculate the excited modes contribution to (A) than the one followed 

in Appendix B. Unlike the proof of [28], however, our calculation may be extended to the case of 
th^L(A) foj. general t in which case matrices {Ani^n'i') in formula (B.l) do not vanish. 


21n2-l / 
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12 Examples 

12.1 Case N = 2 


In the case of two wires, the dimension M of the brane should be 1 for an interesting junction, 
since M = 2 leads to a disconnected junction with S = I and M = 0 gives S = —I, which does 
not conserve the total charge. For M = 1, let k = (a, b) and the compactification radii rf and r|. 
The injectivity of k requires a A b = 1, and the conservation of charge in ensured for a = 5 = 1. 
Forgetting this last requirement for a while, the S'-matrix takes the form 


S = 


1 

rf -I- r^b'^ 


f r\aP' — r^b'^ 
2r|a6 


2rfa6 \ 
r'^b'^ — rfa^J ' 


( 12 . 1 ) 


Two simple but interesting cases arise here. The first one will require the charge conservation 
(a = 5 = 1) but will keep general radii of compactification for each wire, and the second one 
will relax the charge conservation for the equal radii ri = r 2 = r. In the second case, we shall 
consider only the heat transport. 


35 



General radii, charge conserved. Here 


S = 


1 


2r| 


2r1 


O = 


' 1 T / 2 \ “'2 '2 ' 1 ^ 

see (4.29). Note that in the particular case ri = r 2 , 


r? + r'i 


2rir2 


2rir2 


s = o = 


0 1 
1 0 


( 12 . 2 ) 


(12.3) 


corresponding to the fully transmitting junction. For general radii, one obtains from Eqs. (9.13) 
and (9.20) for the mean electric and thermal currents in the non equilibrium stationary state the 
expressions 


^neq('^l (t, 3^)) — ^neq('^2 “ 

^neq(.^l (t, 3j)) = CUneq(.^2 (^5 ^)) 


r^r| V 2 -V 1 
rf + 27r 

(rj + r|)2 V 27r 


(r!Vi+r^,V2) + ^(- 



implying for the electric and thermal conductance the formulae 

eip 1 f-1 ^th 27r rfri /-I l\ 

27r rj + \ 1 “1/ 3/3 (rf + r^)^ \ 1 —1/ 


(12.4) 

(12.5) 


( 12 . 6 ) 


Hence, in mean, (with our convention) the electric current flows through the junction from the 
wire at higher potential to the one at the lower one and, when the potentials are equal, the energy 
current flows from the wire at higher temperature to the one at lower temperature, although the 
latter direction may be reversed by putting the lower-temperature wire in sufficiently high electric 
potential. The large deviation rate function associated to charge only is 




1 r?rl /rj 

2tt {rj + rlY V/3i 



+ 


• 2 2 
1 nrj 

27r r\ -\- ri 


(V2 


Vi)i2, 


(12.7) 


where 12 = — 1 / 2 , see Eq. (10.16). In the special case ri = r 2 = r of a fully transmitting 

junction 

which is compatible^ with Eq. (86) of [11]. The quadratic dependence of on 1 / implies that 

for large time the charge transfers per unit time become Gaussian random variables with mean 
and covariance equal to 


Agi 

t 


Aq2 \ ^ rlrl V2 - Vi 
* / rf +r1 2 tt 


C = 


1 rlrl 

TTt (r? +r2)2 



(12.9) 


To illustrate the latter formulae, we trace in Fig. 4 the dependence of ^(Aqi) and of ^Cn on 

r*! t r 

p = ^ for few values of potential difference and temperatures. 





Figure 4: Mean and covariance for different values of AV = V 2 — Vi, f3i and ^ 2 - 


®Ref. |11] uses a different normalization of the 17(l)-charges so Ui and Vi there are rescaled by ^ relative to 
the ones used here. 
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The large deviation rate function associated to energy only is 


‘V(A) = 



rfrl + rll32Vi)\^ + i/3i/32(Vi - V2){rlVi + rlV2)\ 

2-k + 4ir^r|(/32 - /3i)A + PiP2{t'\ + 

(^l+e-2’r*(/3l+^2))-(Oii)2(e-2’f®/3i+e-2^«/32)-(Oi2)^ ^e-2^*(3l +1^) +e-2’r*(f<2-1^) ) 
^l_g-2,x/3i) (^l_g-2,x/32) 


da; (12.10) 


for A = Ai — A 2 , see Eq. (10.47). For a fully transmitting junction with ri = r 2 = r, the integral 
becomes computable, resulting in the expression 

r^ {PiVl+hV^)\^ + '^Pip2{V?-Vr)X 1__M Hdlll 

Stt (/3i+iA)(/32-iA) ^12V/3i+iA /3i ^2-iA pj ^ ’ 

which agrees with Eq. (90) of [11] taken at ly = 0. Let us look more closely at the analytic 
continuation ***/(—zA) = /(A) of the rate function (12.10) for Ei = E 2 = 0. /(A) is finite for 
—/3i < X < (32 and symmetric around |(/32 — (3i)- Outside that interval, /(A) diverges to + 00 . 
Fig. 5 presents the graph of /(A) and of its Legendre transform 


I{x)= max {Aa;- /(A)} 
Ae]-/3i,/32[ 


( 12 . 12 ) 


for (3i = 1, (32 = 5 and p = — = 1,2,3. The change with increasing p is clearly visible. In 


/(A)=«'/(-iA) Hx) 




Figure 5: Large deviation rate function f(A) and its Legendre transform I(x) for different p = 
r2Mat/3, = l,/3. = 5,l/r=F2 = 0). 

the limit p —> 00 , function /(A) vanishes in the interval ] — /3i,/32[ and stays infinite outside of 
it. The large deviations rate function I{x) is that of the probability distribution of the energy 
change in the first wire per unit time I{x) has linear asymptotes with the slopes 

—/3i and (32 on the left and on the right, respectively, indicating the exponential decay of the 
distribution function of AHi{t)/t arising at long times, with the rate linearly growing with time. 
Fig. 6 zooms on the central region of I{x) around a: = 0 and, for illustrative purpose, presents 




Figure 6: Influence of p on the rate function (zoom), and on the corresponding normalized density 
probability 

the graphs of normalized distribution functions oc exp[—/(a;)]. The influence of p on /'(O) 
and /"(O) representing the long-time mean of AHi{t) and of its variance, both divided by t, is 
depicted in Fig. 7. The mean and the variance per unit time represent, respectively, the mean heat 
current and the thermal noise in the first wire. The increase of p > 1 increases the absolute value 
of the current and decreases the noise. Both exhibit the p ^ tj p symmetry implying that they 
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Mean/t Variance/t 




Figure 7: Influence of p on the mean and the variance of AHi{t) 


are least sensitive to the change of p around p = 1. The influence of the temperature on the rate 
function /(A), its Legendre transform I{x), and on the probability distribution oc exp[—/(cc)] is 
illustrated on Fig. 8 . The asymmetry of the curves increases when the temperature of the second 
wire is lowered below that of the first wire. 



Figure 8 : Large deviations rate function /(A), its Legendre transform I{x), and probability 
density oc exp[—/(a;)] for different P 2 (N = 2, (3i = 1 , r 2 = 2 ri, and Vi =¥2 = 0 ) 


Same radii, charge not conserved. The interest in this case is due to the fact that it 
corresponds to a reflecting and transmitting junction for wires of the same type. Indeed, for 
ri = r 2 = r but a ^ b, 


s = o = 


1 


2a 



for a = - . 

a 


(12.13) 


Since charge is not conserved when a 7 ^ 1, we focus on the energy transport only, and set 
Vi=V 2 = 0. Then 


Wn 





(12.14) 

(12.15) 


and the large deviation function for energy transfer is, with A = Ai — A 2 , 


th 


'/(A) = 


— In 

Jo 


( 1 - 5^2 (l_e-2’rlx>))_e-2irx/32(l_S2^(l_e2irix>.p_|_g-2,rx(/3l+/32) 


dx, (12.16) 


which is illustrated in Fig. 9. The above expression simplifies for a = 1 when S '12 = 1, reducing 
again to the relation ( 12 . 11 ) for fully transmitting junction. 



Figure 9: Large deviations rate function /(A), its Legendre transform I{x), and probability 
density oc exp[—/(a;)] for different a {N = 2, Pi = 1, P 2 = 5, r 2 = ri, and Vi =¥2 = 0) 
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12.2 Case M = 1 


The case M = 1 with the total charge conservation corresponds to the brane B = 
embedded into U{1)^ so that 

U{1) diagonally 

«;= (!,...,!) 

(12.17) 

leading to the S'-matrix 


Sij = PiSij + Ti(l 

(12.18) 

where 

_ 2r] 

Pi ^2 ^ 'Ti 2 

‘^k ' k 2^k ' k 

(12.19) 

are the reflection and transmission coefficients. For the equal radii Vi = 


N-2 2 

Pi — P — jv’ AT 

(12.20) 

leading to a simple nontrivial S'-matrix 


Sij = — p^ij -l- t(1 — 5ij) 

(12.21) 


that was already considered in [36], see also [9]. 


Application to 3 wires We consider the simplest case with the same radii of compactification. 
Here we have p = 1/3 and r = 2/3 and the 5'-matrix is 



( 12 . 22 ) 


The charge and energy currents are 

(^neqiJl) = ^(1/ + I3 - Vi) 

ujneqiKt) =^-^{V2 + Vs- 2Hi)(Hi + V2+Vs) + ^{^ + ^-^) (12.23) 


and similarly on the other wires, cyclicly permuting the indices. For the electric and thermal 
conductance, this gives: 

1 1 \ /-2 1 1 \ 

-1 1 , ‘hg = 1 -2 1 . (12.24) 

1-1/ VI 1 -2/ 

The large deviation function for the FCS of charge transfers is 

®‘./’(j^) = + cycl.)+i^(Vi(2i/i - ^2 - + cych), (12.25) 

where “cycl.” stands for terms obtained by cyclic permutation of the indices. The large deviation 
function for the FCS of energy transfers for 14 = 0 is: 

pOO 

*V(A) = -| In^^dx (12.26) 



where 


N{x, A) = 1 - e-2^^(/5i+/32+/33) p 2 ^g-2,r3:(/3i+/32) ) _ ^2 ^g-2,r3:/3i ^^y^p ^ 

_p .^2 ^g-2'n-a;(/3i+/32) ^g-27ria:(A2-A3) _|_ g-27ria:(Ai-A3)^ _|_ gy^j ^ 

+ Cycl. ) , 

D{x) = (1 - - e-2^^^^)(l - , (12.27) 
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Upon the analytic continuation, *^/(— iA) = /(Ai 2 ,Ai 3 ) for A 12 = Ai — A 2 and A 13 = Ai — A 3 
which is finite only in the region 

— Pi < Ai 2 < P 2 ) —Pi < Ai 3 < Ps , —P 2 < Ai 3 — Ai 2 < Ps ■ (12.28) 

Function / is plotted in Fig. 10 for {Pi, P 2 , P 3 )= (Ijljl) (the equilibrium case) and {Pi,P 2 ,P 3 ) = 
(1,2,3) in the coordinate system with axes at 120° so that the counter-clockwise rotation of the 
graph by 120° corresponds to the cyclic permutation (Ai, A 2 , A 3 ) 1 —>■ (A 3 , Ai, A 2 ). In equilibrium, 
/ is symmetric under such a transformation but out of equilibrium, the above Z 3 symmetry is 
broken to a degree that may be used as a measure of distance from equilibrium. The Legendre 


/(Ai 2 ,Ai 2 ) for (/3 i,/ 32,/33) —(1,1,1) and for (/3i,/32,/33) = (l,2,3) 



-1 -0.5 0 0.5 1 -1.5 -1-0.5 0 0.5 1 1.5 

Ai2 Ai2 

Figure 10: Rate function /(Ai 2 ,Ai 3 ) for 3 wires 

transform of *^/(—iA) is infinite out of the plane xi + X 2 + X 3 = 0 and on that plane, it may be 
regarded as a function 

I{xi 2 ,xi 3 ) = max { 3 ( 2 a:i 2 Ai 2 — 2 : 12 A 13 — 2 : 13 A 12 -I- 2 x 13 A 13 ) ~ /(A 12 , A 13 )} . (12.29) 

Ai 2 )Ai3 

Fig. 11 presents the plot of 1 ( 2 : 12 , X 13 ) for the equilibrium and nonequilibrium choice of tempera¬ 
tures. The level lines of I are equally spaced in various direction far from the origin, indicating the 


T(xi 2 ,xi 2 ) for (/3 i,/ 32,/33) —(1,1,1) and for (/3i,/32,/33) —(1,2,3) 


r^i 


-10 -5 0 5 10 

3^12 

Figure 11: Legendre transform /( 2 ;i 2 ,xi 3 ) for 3 wires 

asymptotic linear increase of the function. The similar breaking of Z 3 symmetry as for / may 
be observed. Finally, Fig. 12 plots for illustration the probability densities oc exp[— 1 ( 2 : 12 , 2 : 13 )]. 



-6 -3 0 3 6 

3:12 


normalized exp[— 7 ( 3 : 12 . 3 : 12 )] for (/3i,/32,/33) = (l,l,l) and for (/3i ,/32,/33) = (l,2,3) 



-6 -3 0 3 6 -6 -3 0 3 6 

3:12 3:12 


Figure 12: Probability density oc exp[—/(a:i 2 , 2 : 13 )] for 3 wires 

Note that most mass of the distribution is in the negative quadrant indicating the heat transfer 
from the hotter 1 ®* and the 2 "'^ wires to the colder 3''“^ one. 
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13 Conclusions 


We have studied a model of a junction of N quantum wires. The Luttinger liquids in the bulk 
of the wires were represented by a toroidal compactification of the A^-component massless free 
bosonic field, with the junction modeled by a simple boundary condition restricting the values of 
the compactified field to a brane B forming a subgroup isomorphic to U{1)^ of the target torus 
U{1)^. The brane B was assumed to be invariant under the diagonal multiplication by phases 
in order to assure the global U{1) invariance and the conservation of the total electric charge. 
We constructed the theory on the classical and the quantum level and showed that the boundary 
condition at the junction leads to a linear relation between the right-moving and the left-moving 
components of the electric currents in different wires. Such a relation describes the scattering by 
the junction of charges carrying the current. The equilibrium state of the system of connected 
wires of length L kept at inverse temperature j3 and in electric potential V was discussed in 
the functional-integral language and in the open-string and closed-string operator formalism, the 
latter being well suited to describe the thermodynamic limit L —>■ oo. We obtained the exact 
solution for the equilibrium current correlation functions both for wires of finite length L and for 
L = oo. In the latter case, the resulting theory provides a special case of the system studied in [32] 
and we adapted from that paper the construction of a stationary nonequilibrium state (NESS) in 
which the wires are kept at different temperatures and electric potentials. Following the lines of 
[10, 11], it was shown that such a state is attained at long times if we prepare disjoint semi-infinite 
wires in equilibrium states at different temperatures and potentials and then connect them by 
the junction and let the dynamics operate. This is a particular realization of the scenario for 
construction of quantum nonequilibrium states proposed in [41]. By considering the constructed 
NESS close to equilibrium, we extracted formulae for the electric and thermal conductance of the 
junction. 

The main result of this paper has been the calculation of the full counting statistics (ECS) for 
charge and heat transfers through the junction and the analysis of its large deviations asymptotics 
at long transfer times in the presence of both transmission and reflection of the conserved charges. 
This was done hrst for charge, then for heat, and finally, jointly for both. We confirmed by an 
exact calculation that the large deviations regime of the charge ECS for a junction of semi-inhnite 
wires may be obtained from a large class of limiting procedures sending to infinity the length of 
the wires as well as the transfer time. The computation of ECS for heat transfers was explicitly 
done only aligning the length of the wires and the time of the transfer, although we developed 
tools for performing the general calculation. We expect that the result for the large deviations 
regime of the ECS for heat transfers obtained in the explicitly treated case applies as well to 
the situation when the length of the wires is sent to infinity faster than the evolution time. The 
expressions obtained for the large deviations rate functions of ECS were compared with the ones 
given by the Levitov-Lesovik formulae. For the charge transfers, we showed that our results for the 
junctions under consideration differ from the Levitov-Lesovik formula for free fermions, although, 
for the vanishing Luttinger couplings, some similarity could be observed in the quadratic part of 
the rate functions that describes the central-limit asymptotics. For the energy transfers through 
the junction, the part of the ECS that was contributed by the excited bosonic Fock-space modes 
appeared to coincide with the bosonic Levitov-Lesovik-type formula for ECS obtained in [28]. 

The simple class of conformal boundary defects considered in the present paper have been 
chosen for illustrative purpose rather than from phenomenological considerations. The latter 
might require introducing a larger family of boundary defects. The simplest extension of the 
class considered here would include conformal boundary defects with displaced branes goB for 
go G or/and the ones with added Wilson lines (w.r.t. a constant gauge held on B). Such 

boundary defects could be dealt with by the same technique, leading to a richer class of S'-matrices 
describing current scattering. More complicated conformal boundary defects would require more 
powerful boundary GET techniques for calculation. For the case of vanishing Luttinger couplings, 
one could use the nonabelian bosonization [48, 6] that comes with a class of conformal boundary 
defects with nonabelian symmetries. Some partial results in this direction have been already 
obtained [23]. The other physically relevant question, not disjoint from the previous one, is the 
stability of the boundary defects in the renormalization group sense. This problem was addressed 
for some simple cases of junctions in [36, 34, 35]. It can be also studied with the boundary GET 
techniques. We postpone a discussion of the above questions to the future research. 
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Appendix A 

Here we shall calculate directly the quantity 


—iu f (0,s) ds 


ujpy(e i ) A-l) 

for one wire of length L > t with the Neumann boundary conditions. In this case, 

(A.2) 

for a 2 n = ^Q^ 2 n and 

t 

[ J'(0,»)d» = + = Ji(/(0,()-/(0,0)), (a.3) 

0 ”5^0 

where the chiral field (p^ is given by (2.21). We shall reorder the exponential of that operator 
writing 
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(A.4) 


Note that the last factor, that may be interpreted as providing the Wick ordering of the left- 
moving vertex operators and e'^*^ (o.o)^ jg ultraviolet singular. We shall replace 

it by its regularized version 
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-Cal 


0 n,<AL = 0 

where A is the ultraviolet cutoff. This leads to the definition: 
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The orthonormal basis of the Fock space Te is given by the vectors 


n 
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(A.IO) 


m—1 


with all but a finite number of pm = 0,1,... equal to zero. Such vectors are eigenvectors of the 

OO 

Hamiltonian with eigenvalues X) and are annihilated by Q^. Hence 
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Since by a straightforward calculation 
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The trace over the zero mode space "Ho spanned by the orthonormal vectors \k) with fc G Z 
such that do|fc) = 2r“^|fc) is 


Multiplying those two traces and the factor of the last line of (A.6) and dividing them by the 
partition function, we obtain the expression 
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By Eq. 16.30.1 of [5], 
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for the first of the Jacobi theta-functions. Hence (A.15) may be rewritten in the form 
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The substitution of the above relation to the (ultraviolet regularized version of) (10.8) results in 
the identity (10.24). 


Appendix B 


Let us consider a quadratic selfadjoint operator 


A = 


N oo 

E E V2 

i.i' —1 n^n' — l 




acting in the Fock space of the vacuum representation of the canonical commutation relations 
(CCR) 

= 0 = ~ ^zz'^n, —n' (®-^) 

built on the normalized vacuum state |0) annihilated by operators a„i. We assume that A„j „/j/ = 
An'i',ni s-iid Bni^n'i' = ,ni- Notc the Commutation relations 
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The exponentiation of the above relation gives 


d-4/4: 

^n: 
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where, in the language of infinite matrices, 
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One of the results of the theory of the Bogoliubov transformations associated to the quadratic 
Hamiltonians A is the formula 

(0|e*-^|0) = det(e*-®P)~^^^ (B.6) 

holding under conditions that guarantee the finiteness of the left and right hand sides which will 
be satished in the cases of interest for us, see e.g. [17]. 

We shall need a generalization of that formula to the expectations of operator e*'^ in certain 
mixed states. Consider such a state Wp corresponding to the density matrix 


N 



E ‘ 
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with tin > 0 diverging sufficiently fast when n —>■ oo and the normalization factor 

N oo 
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We would like to calculate the expectation value a;p(e®'^). To this end, we may use the Araki- 
Woods representation [3] of the CCR acting in the double Fock space F ® F and given by the 
formula 

I_V + bni iahiY + VihHani\ iR Ql 

[cinij \dm/ I yi + bni + \/b^ V ' 


where 


^ , d-ni 


= ani® I, = I ® an 
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The Araki-Woods representation has the property that the matrix element on the vacuum |0,0) = 
|0) ® |0) of a CCR observables taken in that representation reproduces their Wp expectations. In 
particular 

Wp(e'-^) = tr^(^pe'-^^ = (0,0|e*-^|0,0) , (B.ll) 


where 
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we infer from the previous result (B.6) that 
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The above calculation can be applied to the operator A corresponding to the excited modes 
part of ^ \ for any time t, but we shall limit ourselves to the simpler instance 

I 

when t = 2L. In that case 

00 
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see (10.41), with the identification 
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For the density matrix, 

we shall take 
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so that the state ojp coincides with ojq on the algebra generated by the excited modes. The 
normalization 

Z0 = (B.20) 

i,n>0 

In this case = (e^”^* — 1)“^. It follows that 
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A little of straightforward algebra shows that 
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Hence, denoting by &„ the diagonal N x N matrix with entries bni, we obtain 
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On the other hand, 


det (e*'®"’") 

We infer then from (B.15) that 
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What we have to compute, however, is the contribution of the excited modes to (10.42) which is 
equal to 

E a(-2n)ia(2„)i .\ 

i n>0 ' 


itr(e 


(B.26) 


This is clearly obtained by multiplying (B.25) by shifting /3 to /3 + iX in the result, and 
re-dividing it by Z^, which gives 
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i.e. the result (10.45). Note that the last expression is invariant under the change of A to —A+i/3. 
Indeed, 
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Appendix C 

We compute here the scaling limit 

lim 6»2(^(6»-izy)| (C.l) 

where the rate function (j){6~^u) is given by relation (11.2). The logarithms of determinants on 
the right hand side of (11.2) will be computed by expanding: 
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Al± = /±(e)|^_,_,^(eT‘^~'‘^Ste±i«-^-S-/) = /±(e)|^ ,_,^(±id-ii? + 0-2i? + O(d-3)) (C.3) 
with 
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Only the first two terms of that expansion will contribute to the scaling limit (C.l). The integra¬ 
tion over e is reduced to the terms 
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Thus 
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Using formulae (C.4), we finally get: 
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The other non-vanishing terms in the scaling limit (C.l) come from 



00 

J tr(A±)de = - E A') , 

0 *■*' 

(C.9) 
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The distinction between the contributions for different signs disappears at this order 
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On the other hand, for i ^ i', 
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The addition of contributions (C.8) and (C.15) results in the identity (11.3). 
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